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Abstract 


Background: Arrhythmogenic Cardiomyopathy (ACM) is a potentially 
lethal form of cardiomyopathy, with a high incidence of mortality. It 
is prominently characterised by abnormal ventricular anatomy, fibrofatty 
replacement of the myocardium and the His-Purkinje system, as well as 
electrophysiological remodelling. It remains unknown how the broad spec- 
trum of ACM-related co-morbidities alter the 12 lead electrogram (ECG). 
Ventricular activation abnormalities result in QRS complex and terminal 
activation duration prologation, as well as abnormal R-wave progression 
and low amplitude signals are diagnosis criteria of ACM on the 12-lead 
ECG. The pathological remodelling leading to the above electrocardio- 
graphic findings can be of electrical, geometrical or structural nature or 
a combination of substrates. Diagnosis of this clinical entity is challeng- 
ing and relies on a combination of ECG, non-invasive imaging approaches 
to access myocardial muscle structrure and geometry as well as genetic 
testing. Despite the multi-modality diagnosis approach, misdiagnosis is 
common and some patients are not diagnosed until post-mortem exami- 
nation is performed. Improvements in diagnosis are much needed and one 
challenge is understanding the link between ECG manifestations and un- 
derlying disease burden. Furthermore, no ACM specific treatments exist 
to date and life quality of sufferers is greatly reduced as cardiac defib- 
rillators and guidelines to avoid physical activity are the treament norm. 
Obtaining an accurate diagnosis non-invasively can contribute to better 
risk stratification and patient management, ultimately improving the life 


quality and safety of patients. 


Aims: It is known that several mechanisms can lead to slow electri- 
cal conduction in the ventricles of ACM patients. Plausible single cell 
electrophysiological mechanisms, such as co-modulation of Iyq-I«1 ionic 
channels, leading to tissue level electrical conduction velocity modulation 
in ACM are explored. This thesis aims to uncover which mechanisms of 


conduction slowing prevail in patients with ACM by integrating informa- 
tion from imaging data and performing personalised activation sequence 
simulations of patients diagnosed with ACM. The scope is to also relate 
the different conduction abnormalities to signature traits on the ECG, so 


as to help refine diagnosis and disease burden non-invasively. 


Methods: The effect of ACM specific ionic channel joint reduction of 
conduction is explored via a study relating co-modulation in key ionic 
channels (Iya-L1) [1], known to be remodelled in ACM variants|2| and 
the tissue level electrical conduction velocity in 2D tissue simulations. 
Personalised models of electrical conduction are created from Magnetic 
Resonance Imaging (MRI) data of a cohort of ACM patients with 15 
definite, 3 borderline and 2 plausible patients according to the 2010 inter- 
national diagnosis guidelines [3]. Patient ventricular and torso anatomies 
are extracted from the MRI images and regional fibrotic scar informa- 
tion is extracted from contrast agent enhanced MRI imaging. Different 
types of remodelling were introduced into graph-based models of sinus 
rhythm activation, namely localised fibrosis-mediated conduction slowing 
in the myocardial muscle or in the Purkinje-myocardial fast conduction 
system, myocardial tissue uncoupling leading to global ventricular tissue 
conduction slowing, right ventricular (RV) muscle disappearance leading 
to slowed conduction localised in the RV, RV muscle volume dilation and 
RV wall thinning. Insights from slab tissue simulations were used to es- 
timate the order of magnitude of conduction velocity modulation due to 
ionic remodelling, clinical data (contrast enhanced MRI) for scar based 
conduction anomalies and from other computational studies quantifying 
conduction velocity modulation in the ACM myocardium [4] were used to 
find physiological parameter intervals for the computational study. The 
effects of the above types of remodelling on features of the ECG relating 


to activation dynamics were quantified. 


Results: Slab tissue simulations revealed that ACM-specific ionic re- 
modelling of the Iyg-I«¢, macromolecular complex can result in up to 
45% bulk myocardial tissue CV modulation. A computational pipeline 
was adapted to access the effect of ACM ventricular substrates on con- 
duction properties of the ventricles. The personalised data-driven models 
of two ACM patients with varying disease burden were used to success- 


fully explain how the underlying whole ventricular remodelling impacted 


the clinical patients’ ECG recordings. Computer simulations demonstrate 
that, in ACM, mild QRS prolongation (+100 ms), in absence of other QRS 
abnormalities, is explained by regional myocardial conduction slowing in 
the areas of MRI-identified fibrosis. However, MRI-segmented myocar- 
dial fibrosis and RV dilation alone proved insufficient to recapitulate more 
severe QRS abnormalities, such as severely prolonged mean QRS dura- 
tion (>120 ms) across precordial leads, mean terminal activation dura- 
tion across V1-V4 (>55 ms), and low QRS mean R-wave amplitude (<2 
mV) in precordial leads, as well as poor R-wave progression across precor- 
dial leads. A 70% decrease in Purkinje-endocardial conduction and 30% 
RV wall thinning were additionally required to explain these more severe 


abnormalities. 
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Chapter 1 


Introduction 


In this chapter, the relevant biological concepts underpinning and motivating this 
work, as well as ACM, the disease of interest, are presented. The evolution of the 
clinical disease in time is followed by the basic functioning of the heart and the role 
of ventricular electrophysiology in determining the action potential and maintaining 
electrical conduction. Then, the ways in which ACM changes ventricular electro- 
physiology are explained with reference to known genetic mutations leading to the 
experimentally observed disease manifestations. Hence, a link between extracellular 
level mechanisms and global structural alterations to the cardiac muscle properties 
is presented in the context of known genetic mutations of ACM. Finally, the causes 
of conduction slowing, leading to abnormal ventricular electrical activation found in 


this disease are summarised. 


1.1. Arrhythmogenic Cardiomyopathy - A novel clin- 
ical entity 


1.1.1 Historical Perspective of the clinical evolution of Ar- 
rhythmogenic Cardiomyopathy 


Medical practice has continuously updated the definition of the clinical entity of 
interest, with its currently accepted name being Arrhythmogenic Cardiomyopathy. 
Back in the 1980s [5], the first reports by clinical doctors were of a muscle disease 
affecting the right ventricle, ultimately leading to muscle damage and lethal ventric- 
ular arrhythmias. Since then, numerous genetic and developmental factors have been 
outlined and the clinical entity has evolved to encompass both ventricles. Strict di- 
agnosis criteria have been proposed in 1994, which were then updated in 2010 and 


2019 respectively. The progression has been towards finding patient subpopulations 


and classifying patients in risk categories, based on an extensive, non-invasive clin- 
ical examination. Clinical attention has shifted from the right ventricle only to a 
bi-ventricular disease, whereby either of the two ventricles can be remodelled and 
an arrhythmic subsrate can be a result of muscle degradation and subsequent fibro- 
fatty replacement. The name has evolved from Right Ventricular Dysplasia [5] to 
Arrhythmogenic Right Ventricular Cardiomyopathy (ARVC) [6] to the current name 
Arrhythmogenic Cardiomyopathy (|7]. 


1.1.2 ARVC - A right ventricular dominant clinical entity 


Sometimes the literature mentions the name Arrhythmogenic Right Ventricular Car- 
diomyopathy, when referring to a heart disorder that primarily manifests itself in the 
right ventricle in humans. It affects the muscle tissue in the heart, also known as 
myocardial tissue. It is characterised by a progressive decoupling of the myocyte net- 
work. As a result, micro disarrangements, fatty replacement and fibrosis characterise 
the muscle tissue and create a pro-arrhythmic layer. Thus, ARVC can lead to ventric- 
ular tachycardias, ventricular arrhythmias and/or sudden cardiac death [8] [9] [10]. 
ARVC is a silent killer as micro-scale structural and localised electrophysiological ab- 
normalities are often not visible in clinical examinations and these can be sufficient 
to lead to lethal arrhythmias or sudden cardiac death (SCD) in certain patients [11] 
[12] [5] [13]. In particular, young patients and athletes are at high risk [13] [14]. It 
is worth noting that ARVC is known as Arrhythmic Right Ventricular Dysplasia in 
the US. ACM prevalence is between 1:2000 and 1:5000 in the general public in the 
US [15] and affects men more frequently than women, with a ratio of 2.4:1 [16]. The 
disease has no specific racial or geograghical predisposition 

The disease is diagnosed universally according to the 2010 task force criteria [3], 
which build and refine the 1994 criteria [6]. A comparison of the criteria as they stood 
in 1994 and their update in 2010 [3] is included in table 1.1 and 1.2. The most typ- 
ical clinical diagnosis has 2 main scenarios: (1) electro-cardiographic abnormalities, 
typically a ventricular tachycardia with left bundle branch morphology and T-wave 
inversion in leads V1-V3 (2) functional and structural abnormalities mostly of the 
RV, such as wall thinning, regional wall motion alterations and global dilation. The 
disease is clinically heterogeneous, occurring with inter-famillial and intra-famillial 
variability. The left ventricle is frequently involved in the disease and hence ACM is 
a more useful clinical name for this disease [17]. This is an umbrella term, used to 
describe both arrhythmogenic left ventricular cardiomyopathy and right sided man- 


ifestations of disease. This new clinical entity is currently under investigation, with 


leading clinicians currently developing novel scoring systems for diagnosis based on 
MRI, ECG, histology and genetics [18] [19]. In this thesis, the nomenclature ACM 
will be used throughout except for in cases when the observations are solely valid for 


the right ventricular manifestation of disease. 


1.2. Basic Function of the Heart 


The heart is a mechanical pump and is used to pump blood around the circulation 
system [20]. This is achieved by the heart muscle doing active work on the blood in 
the chambers of the heart. The heart is separated in two atrial chambers used for 
filling and two ventricular chambers used for pumping the blood into the two major 
arteries. In order to achieve muscle contraction, a wave of electricity propagates along 
the heart muscle. Physiological propagation occurs when the stimulus is dictated 
by the right atrial sinus node pacemaker and by the atrio-ventricular node for the 
ventricles. The atrial node cells have automaticity, meaning they are self-excitng [21]. 
The propagation then descends the Purkinje network of fibers, which is connected to 
the endocardium of the ventricles and from there proceeds to reach all parts of the 
ventricular muscle. This process is illustrated in Figure 1.3. On the mechanical 
side, atrial contraction pumps blood into the ventricles, which in turn apply pressure 
to pump it into the Pulmonary artery and the Aorta in a coordinated and timely 
manner. Once the contraction of the whole heart has been triggered by the electrical 
stimulus, the heart enters in a relaxation stage when the atrial chambers refill as the 
four chambers regain their end of contraction volume, also known as the diastolic 
volume. 

The propagation of the electrical wave through the heart muscle, also known as 
the myocardium, is maintained by a flow of ions between the cardiac muscle cells via 
a coupling structure known as the gap junction [23]. Heart muscle cells are arranged 
in fibres [24]. There is intrinsic anisotropy in heart muscle fiber propagation speed. 
This arises as electrical conduction speed along the fiber direction is higher than speed 
along the transverse or interfiber directions. As the electrical impulse travels from 
cell to cell, it causes a cyclic change of ionic concentrations, which give rise to the 
action potential, a potential difference of the cell with respect to the extracellular 
space [25]. 


Original Task Force Criteria Revised Task Force Criteria 
I. Global or regional dysfunction and 
structural alterations” 


Major 
By 2D echo: 
* Severe dilatation and reduction of RV * Regional RV akinesia, dyskinesia, or aneurysm 
ejection fraction with no (or only mild) LV * and 1 of the following (end diastole): 
ampeirment, — PLAX RVOT =32 mm (corrected for body size [PLAX/BSA] =19 mm/m2) 
* Localized RV aneurysms (akinetic or — PSAX RVOT =36 mm (corrected for body size [PSAX/BSA] =21 mm/m2) 
dyskinetic areas with diastolic bulging) 
+ Severe segmental dilatation of the RV — or fractional area change <33% 
By MRI: 
« Regional RV akinesia or dyskinesia or dyssynchronous RV contraction 
* and 1 of the following: 
— Ratio of RV end-diastolic volume to BSA =110 mL/m? (male) or =100 
mL/m? (female) 
— or RV ejection fraction <40% 
By RV angiography: 
* Regional RV akinesia, dyskinesia, or aneurysm 
Minor 
By 2D echo: 
* Mild global RV dilatation and/or ejection + Regional RV akinesia or dyskinesia 
fraction reduction with normal LV * and 1 of the following (end diastole): 
— PLAX RVOT =29 to <32 mm (corrected for body size [PLAX/BSA] =16 to 
* Mild segmental dilatation of the RV <19 mm/m2) 
* Regional RV hypokinesia — PSAX RVOT 232 to <36 mm (corrected for body size [PSAX/BSA] =18 to 


<21 mm/m?) 
— or fractional area change >33% to <40% 
By MRI: 
* Regional RV akinesia or dyskinesia or dyssynchronous RV contraction 
+ and 1 of the following: 
— Ratio of RV end-diastolic volume to BSA =100 to <110 mL/m? (male) or 
=90 to <100 mL/m?2 (female) 
— or RV ejection fraction >40% to <45% 
II. Tissue characterization of wall 


Major 
* Fibrofatty replacement of myocardium on * Residual myocytes <60% by morphometric analysis (or <50% if 
endomyocardial biopsy estimated), with fibrous replacement of the RV free wall myocardium in =1 
sample, with or without fatty replacement of tissue on endomyocardial 
iopsy 
Minor 


+ Residual myocytes 60% to 75% by morphometric analysis (or 50% to 65% 
if estimated), with fibrous replacement of the RV free wall myocardium in 
21 sample, with or without fatty replacement of tissue on endomyocardial 
iopsy 

III. Repolarization abnormalities 


Major 
* Inverted T waves in right precordial leads (V1, V2, and V3) or beyond in 
individuals >14 years of age (in the absence of complete right bundle- 
ranch block QRS =120 ms) 
Minor 
+ Inverted T waves in right precordial leads + Inverted T waves in leads V1 and V2 in individuals >14 years of age (in 
(V2 and V3) (people age >12 years, in the absence of complete right bundle-branch block) or in V4, Vs, or Vg 


absence of right bundle-branch block) 


* Inverted T waves in leads V1, V2, V3, and Vy in individuals >14 years of 
age in the presence of complete right bundle-branch block 
(Continued) 


Figure 1.1: The 1994 and 2010 task force criteria for the diagnosis of ACM. The 
adaptation incorporates the newest imaging techniques in the process of diagnostics. 


Original Task Force Criteria Revised Task Force Criteria 
IV. Depolarization/conduction abnormalities 
Major 
¢ Epsilon waves or localized prolongation (>110 * Epsilon wave (reproducible low-amplitude signals between end of 
ms) of the QRS complex in right precordial leads QRS complex to onset of the T wave) in the right precordial leads (V, 


(V1 to V3) to V3) 
Minor 
* Late potentials (SAECG) ¢ Late potentials by SAECG in =1 of 3 parameters in the absence of a 


QRS duration of =110 ms on the standard ECG 

¢ Filtered QRS duration ({QRS) =114 ms 

¢ Duration of terminal QRS <40 pV (low-amplitude signal duration) 
238 ms 


¢ Root-mean-square voltage of terminal 40 ms <20 pV 
¢ Terminal activation duration of QRS =55 ms measured from the 


nadir of the S wave to the end of the QRS, including R’, in Vj, V2, or 
V3, in the absence of complete right bundle-branch block 
V. Arrhythmias 
Major 
¢ Nonsustained or sustained ventricular tachycardia of left bundle- 
branch morphology with superior axis (negative or indeterminate QRS 
in leads II, III, and aVF and positive in lead aVL) 
Minor 
¢ Left bundle-branch block-type ventricular ¢ Nonsustained or sustained ventricular tachycardia of RV outflow 
tachycardia (sustained and nonsustained) (ECG, configuration, left bundle-branch block morphology with inferior axis 
Holter, exercise) (positive QRS in leads II, III, and aVF and negative in lead aVL) or of 


unknown axis 


¢ Frequent ventricular extrasystoles (>1000 per »* >500 ventricular extrasystoles per 24 hours (Holter) 
24 hours) (Holter) 


VI. Family history 


Major 
* Familial disease confirmed at necropsy or « ARVC/D confirmed in a first-degree relative who meets current Task 
surgery Force criteria 
« ARVC/D confirmed pathologically at autopsy or surgery in a first- 
degree relative 
* Identification of a pathogenic mutation! categorized as associated or 
probably associated with ARVC/D in the patient under evaluation 
Minor 
¢ Family history of premature sudden death (<35 * History of ARVC/D in a first-degree relative in whom it is not 
years of age) due to suspected ARVC/D possible or practical to determine whether the family member meets 
current Task Force criteria 
¢ Familial history (clinical diagnosis based on ¢ Premature sudden death (<35 years of age) due to suspected 
present criteria) ARVC/D in a first-degree relative 


« ARVC/D confirmed pathologically or by current Task Force Criteria 
in second-degree relative 
PLAX indicates parasternal long-axis view; RVOT, RV outflow tract; BSA, body surface area; PSAX, parasternal short-axis 
view; aVF, augmented voltage unipolar left foot lead; and aVL, augmented voltage unipolar left arm lead. 
Diagnostic terminology for original criteria: This diagnosis is fulfilled by the presence of 2 major, or 1 major plus 2 minor 
criteria or 4 minor criteria from different groups. Diagnostic terminology for revised criteria: definite diagnosis: 2 major 
or 1 major and 2 minor criteria or 4 minor from different categories; borderline: 1 major and 1 minor or 3 minor criteria 
from different categories; possible: 1 major or 2 minor criteria from different categories. 
*Hypokinesis is not included in this or subsequent definitions of RV regional wall motion abnormalities for the proposed 
modified criteria. 
tA pathogenic mutation is a DNA alteration associated with ARVC/D that alters or is expected to alter the encoded 
protein, is unobserved or rare in a large non-ARVC/D control population, and either alters or is predicted to alter the 
structure or function of the protein or has demonstrated linkage to the disease phenotype in a conclusive pedigree. 


Figure 1.2: The 1994 and 2010 task force criteria for the diagnosis of ACM. The 
adaptation incorporates the newest imaging techniques in the process of diagnostics. 
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Figure 1.3: The principal components generating and carrying an electrical impulse 
throughtout the heart are shown, along with the four chambers of the heart [22]. 


1.3. Single cell-level mechanisms of electrical conduc- 
tion 


1.3.1 Cardiac Action Potential 


The electrical wave of excitation illicits a cyclic change in the ionic concentrations of 
key ions in the cardiac muscle cell or myocyte. More specifically in the ventricles, 
this generates a cell action potential (AP) with a characteristic curve, as outlined in 
Figure 1.4. The ventricular myocytes undergo a distinct 4-phase cycle: 


e Phase 0 Upstroke or otherwise known as rapid depolarisation. The membrane 
voltage shifts to a positive value of +40mV from close to -90mV in resting 
state as positive ions, Nat and Ca?* enter the cell through voltage-gated ion 


channels. 


e Phase 1 Early rapid repolarisation through the opening of potassium repolar- 


ising channels. 


e Phase 2 A plateau phase, dependant on balance of inward Ca?* channels and 


outward repolarising channels. 


e Phase 3 Final repolarisation that brings membrane potential to its resting, 


negative value. 
e Phase 4 Resting potential of approximately -90 mV. 


The electrical wave first leads to depolarisation comprises of ion dynamics entering 
the cell is governed by a cascade of channels, initiated by the rapid outward Nat 
channel. 

The Na* Channel 

As the neighboring cell depolarises, the potential of the current cell increases to 
above its resting value of -90 mV as an inward fast Nat, or NaV1.5 channel opens 
and thus increases the cell’s potential further, until the threshold value of -70 mV is 
reached and a self sustaining inward fast Na* current onsets. The large Nat current 
rapidly brings the trans-membrane potential (TMP) to about 0 mV for a transient 
period of time called the overshoot. L-type Ca?*+ open at -40 mV and initiate the 
calcium cycle described in Section 1.3.1. The fast Nat channel opening is responsible 
for the initial upstroke in the action potential, as shown in Figure 1.4. This channel 
also has a role in conduction velocity modulation, as discussed later in this chapter 


and in the simulations in Chapter 2. 
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Figure 1.4: A cardiac muscle cell’s major ionic currents are responsible for the ac- 
tion potential shape [26], including the delayed rectifier potassium current (1), the 
sodium channel (Iy,), the L-type calcium current (Icqz), the transient outward potas- 
sium channel (J;), the ultra rapid potassium channel (Ixu,), the delayed potassium 
current Ix,, the sodium-calcium exchanger (NCX), the funny current (J;) and the 
muscarinic potassium channel (I«cn). 


The characteristic AP curve varies across the myocardial wall. There is disagree- 
ment in the literature about the exact pattern governing myocardial AP variability 
but gradients in the duration of phases 1, 2 and 3 have been observed experimentally 
[27][28], whilst later wstudies claim that no significant repolarisation gradients exist 
within a healthy myocardium [29]. The duration of the cardiac muscle cell cycle plays 
a vital role in arrhythmia initiation and maintenance, especially in the context of AP 
heterogeneity across one or several of the heart’s chambers. 

The Calcium Cycle 

When a neighbouring cell is excited, the associated potential gradient causes an 
influx of Ca?* into the T-tubules of the currently resting cell. The voltage gated 
Ca?* channel opens in the cellular membrane of the tubules and calcium enters the 


cell. Ca?* induced Ca?* release occurs via the ryanodine receptor channels (RYR2) 


[30]. These extra calcium ions rapidly leave the sarcoplasmic reticulum (SR), where 
they are stored internally to the cell. The abundance of calcium in the cell binds to 
the troponin and initiates muscle contraction through shorterning of the sarcomeres. 
Muscle relaxation occurs as Ca?*+ unbinds from troponin. Ca? is transferred back 
into the SR via a Ca?*+ ATPase, at the expense of ATP hydrolysis. This hydrolysis 
transforms chemical energy, stored as phosphoanhydride chemical composition, into 
useful mechanical work. Calcium is stored in the SR until the next contraction. 
To further reduce the concentration of Ca?*, the remaining intra-cellular Ca?* is 
exchanged with Nat via the sodium-calcium exchanger. As 3 Na* ions enter for every 
expelled Ca?", this contributes to a buildup of a net positive Ca?* inward gradient of 
the cell with respect to the extracellular space. Finally, the Na* gradient is controlled 
via the sodium-potassium pump. Na*/ K* pump is a solute pump that pumps 
Nat out and K* in, against their concentration gradients. The intracellular Nat 
concentration needs to be smaller than the extracellular Nat concentration during 
resting state, as the inward Nat current plays an important role in initiation of muscle 
contraction [31]. A schematic of the processes above is shown in Figure 1.5. 
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Figure 1.5: The calcium ion cycle is central to the governance of muscle contraction 
and relaxation. There is an interplay between the difference in the relative intra- 
cellular and extracellular Ca?* concentrations. This illustration is reproduced with 
permission from Chapter 9 of the Human Physiology : from cells to systems [32]. 


T-Type Calcium and L-Type Calcium channels 

The L-type calcium channel has a long lasting activation in a cardiac cycle and is 
part of the high voltage activated family of channels. It is the link between excitation 
and contraction of myocytes. The inward flow of ions into the cell is responsible 
for triggering the RYR2 receptor on the SR [33]. Phosphorylation of this channel 
increases its permeability and consequently the myocyte contractility. 

The T-type Ca?* opening is responsible for initiating L-type calcium channels 
opening and is transient in nature [32]. It opens during membrane depolarisation, 
mediating Ca?* ion influx via the cell membrane. The T-type channel is activated in 
the same voltage range as the voltage-gated Na* channel, at about -55 mV. The large 
potential gradient induces a considerable flow of ions into the cell. This channel is 


regulated by dopamine and other neurotransmitters, which block T-type currents [34]. 


The Ryanodine Receptor 2 channel (RYR2) 

RYR2 is a protein involved in the process of Ca?* induced Ca?* release in my- 
ocytes. It opens during systole and is located on the SR. It is activated by the binding 
process of Ca?* to the SR wall, with Ca?* receptors governing the channel open and 
closed states [35]. Their activation is Ca** gated and their opening rate and duration 
is dependent on the intracellular Ca?* concentration. 

The Na*/K*ATPase Pump 

The resting state potential of the cell is maintained at -90 mV and is controlled by 
potassium ion dynamics. The Nat/K*ATPase pump moves 2 K* in and 3 Nat out 
at each cycle [36]. It thus contributes to generating a negative resting cell potential 
by creating an effective outward current, whilst building up a positive K* gradient 
between the cell and its surroundings. In the short term, more Kt will exit the cell 
via this gradient than will enter. A negative trans-membrane voltage grows. This 
process also builds a negative Na* gradient across the cell membrane. 

Voltage Gated Kt Channels 

Voltage-gated K* channels are trans-membrane channels that respond to voltage 
changes across the myocyte membrane. 

Transient Outward I,, 

It contributes to early repolarisation phase 1 of the AP by moving intracellular 
K* to the extracellular space. It has rapid activation and deactivation [37]. It is 
activated as a result of the fast increase in membrane potential in phase 0. 

Ultra Rapid Potassium Channel Ix,,, 
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This ultra rapid outward channel opens in the plateau phase 2 of the AP [38]. It 
plays an important role in determining the contractility of the atrium but not that 
of the ventricles. 

Delayed Rectifier I, 

The delayed potassium rectifier helps to repolarise the cell by expelling K* ions 
from the cell in phase 3. It is activated by the influx of Na* during upstroke. It 
contributes to defining the shape of the plateau in the AP. 

Inwardly rectifying potassium channel Ix, 

This inward rectifying current stabilises the resting membrane potential and switches 
on during the initial depolarisation and final repolarisation of the AP. Inwardly rec- 
tifying Kt channels conduct ions into the cell at potentials more negative than the 
Nernst potential, which is the equilibrium potential of this channel, Ex. At equilib- 
rium, the diffusion gradient across the cell membrane exactly opposes the electrical 
charge density within the cell and thus no charge flows through the channel. By doing 
so, they contribute to stabilising the resting cell potential of the cell. These inward 
rectifiers also conduct less outward current, when V > Ex, than traditional voltage 
gated channels. 

Na*/Ca?*+ Exchanger 

The sodium-calcium exchanger, or NCX, is an exchanger protein which removes 
calcium from the myocyte. It uses the energy stored in the eletrochemical gradient of 
Nat, which is built up at the beginning of the AP [39]. This protein takes advantage 
of the high extracellular concentration of Nat and allows Na* to diffuse into the cell, 
in exchange for C'a?* ions out, in the ratio 3:1. Thus, a net increase of positive charge 
in the cell occurs via this mechanism. 

ATP-Sensitive K* Channel 

Karp is a potassium channel gated by intracellular nucleotides, ATP and ADP. 
This channel is found in the plasma membrane or within sub nuclear membranes. 
Absence of Karp significantly impairs the myocyte’s ability to distribute Ca?*, de- 
creasing sensitivity to sympathetic stimulation [39]. 

The slow repolarising / ks channel 

The /K‘s channel is a major voltage gated repolarising current in the ventricles, 
expelling [A ]* ions from the myocyte in phases 2 and 3 of the ventricles [40]. It has 
slow opening kinetics and has been associated with heterogeneous expression across 
the myocardial wall. 

Sarco/endoplasmic reticulum Ca** -ATPase (SERCA) Pump 


it 


The pump uses ATP (chemical stored energy) to transfer Ca** ions back into 
the SR after the contraction has terminated and Ca?* has unbinded from troponin 
in the cell. It regulates cytosolic Ca?* concentration, SR Ca?* load, and the rate 
of contraction and relaxation of the heart [41]. It is controlled by the regulatory 
protein phospholamban in the ventricles and is potentiated under beta-adrenergic 
stimulation. 

Beta adrenergic receptors 

An important alteration of calcium handling is mediated by catecholamines (adrenaline 
and neuradrenalin) belonging to the sympathetic nervous system, which controls the 
fight or flight response. The catecholamines bind to beta-receptors in the heart. Acti- 
vation of these receptors increases heart rate, relaxation speed and contractility. The 
increase in relaxation speed and contractility are thought to be the effect of cAMP de- 
pendent protein kinase (PKA)-mediated phosphorylation of proteins involved in Ca?* 
handling [39]. Three important calcium handling proteins that are phosphorylated 


are: 


e The L-type Ca?* channel, involved in letting calcium into the cardiac cell at 
the beginning of the contraction. The current is enhanced by increasing the 


mean channel open time [42]. 


e RYR2, which controls contractility by gating the flow of sarcoplasmic calcium 


to the cytosol. When phosphorylated, the SR Ca?* release is amplified. 


e The phospholamban (PLB), which normally inhibits the SERCA pump. When 
phosphorylated, this ability is lost. This enables the excess intracellular Ca?+ 
to be transported back to the SR in diastole. 


1.4 Effect of ACM gene mutations on single cell elec- 
trophysiology 


Overwhelming evidence suggests that electrical abnormalities precede clinically per- 
ceivable structural abnormalities in ACM [43] [44] [45]. Examining the genetic basis 
of ACM can help to shed light on the pathogenesis, in the context of electrical disease 
progression. 

Calcium handling is shown to be disrupted in patients with a PKP2 mutation in 
human, cellular reprogramming in vitro studies [46] [47] as shown in Figure 1.6. The 


cells are engineered using patient fibroblasts and the iPSC technique, which is used to 
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induce pluripotent stem cells from fibroblasts or embryonic cells [48]. Patient specific 
cells harbouring a PKP2 mutation with various pathogenic conditions were created in 
an vitro adult like feeding mechanism, with attenuated PPARG expression to produce 
fatty infiltration conditions. The cells were sustained in this environment, which 
simulates a system of adult cardiomyocytes outside of the human body [46] in ACM 
type extracellular conditions. Using this setup, it was shown that the intracellular 
Ca** concentration remains abnormally high in PKP2 mutated cells, as in Figure 
1.6. 
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Figure 1.6: Ca?* imaging in a PKP2 mutated human cell line of ACM shows Ca?* 
cellular overload. (a) The pacing protocol. Myocytes are paced 10 times at 1Hz. (b) 
The intracellular Ca?* decay is shown for 2 different solutions, one having a nor- 
mal 1.8 mM concentration of extracellular calcium and one with abnormally high 
concentration of extracellular Ca?*. Both measurements occur at normal lipogenesis 
and thus mimic normal PPARG functioning. (c) Healthy myocytes (1) and PKP2 
dysfunctioning myocytes (2) were placed in an environment mimicking PPARG dys- 
function and high lipogenesis. The Ca?* decay constant for PKP2 modified cells is 
significantly larger than in healthy myocytes, suggesting a disrupted calcium handling 
is prominent in ACM [46]. 


In order to explain the calcium mishandling further, it is useful to look at levels of 
Ca?* handling proteins. These are shown to be different from healthy myocyte levels 
and are compared in Figure 1.8. 

A cell culture of iPSC-derived cardiomyocytes can be made with an ACM type 
PKP2 gene mutation. In a particular study, fibroblasts from a 30-year old male were 
used to populate the cell cultures and the cells were voltage clamped. The L-type 
channel’s current response to a fixed voltage is measured and shown in Figure 1.7. 


An in vitro analysis of protein expression shows several signalling proteins of ionic 
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channels are disrupted in this mutation, which makes up to 2% of all ACM mutations 
worldwide. The protein expressions change in a fatty environment and can be seen 
in Figure 1.8. In particular, the protein expression of the L-type calcium channel is 
overexpressed. 

A PKP2 knock out mouse model showed a prolongation of Ca?* transient duration 
at 30% recovery and at 50% recovery from peak Ca?* by 25% and 10.5% respectively 
[49]. The APD90 was also increased by 25% in the PKP2 mouse model and reflected 
the increase in Ca?* transient duration. Major Ca?+ handling dysruption is caused 


by the lowered Jcaz peak current associated with an increase in its inactivation time. 
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Figure 1.7: A representative current vs fixed voltage clamp response of Iggz. The 
test voltage pulses in the range from -40 to +50mV. The cell response is that of a 
genetically engineered ACM mutation in the PKP2 gene, a prominent ACM mutation. 
This response is deemed similar to a healthy ventricular myocyte [47]. 
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Figure 1.8: The same preparation protocol is used as in Figure 1.6. (a) In a normal 
PPARG environment, the RYR2 and Cav1.2 handling proteins are over expressed in 
PKP2 mutated cells, where as NCX1 and SERCA2a are under expressed. Note that 
these are merely the proteins handling the channels and not the channel currents 
themselves. (b) In a fatty environment, the protein levels follow a similar trend to 
the one in (a) with a reversal in tendency for the RYR2 protein. Furthermore, the 
change in channel protein expression between healthy and mutated myocytes grows 
in this environment. (* implies p <0.05) 
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Ina is targeted by ACM genetic mutations. In particular, a mutation of the DSG2 
gene in transgenic mice is reported to decrease Nat current density. This results in 
a reduction in the AP upstroke velocity, slower conduction and increased arrhythmia 
susceptibility at disease stages prior to fibrotic replacement and necrosis [50]. The 
corresponding mouse model ionic behavior is illustrated in Figure 1.9. A similar 
behaviour is recorded in genetically engineered human cells, plagued by a SCN5A 
mutation and presented in Table 1.10. 

To further quantify the degree of sodium channel reduction in ACM patients 
and electrophysiological consequences, AP biomarkers from a PKP-2 ACM human 
mutation, studied in iPSC in vitro human cells are included in Figure 1.1. 


Cell Type APA(mV) | MDP(mV) | cAPD90 | n_| 
H9 hESC-CMs 16.942.2 | 60.8+25 | -437418 | 201.8439 


ARVD iPSC-CMs 49.4440 | -37.5+3.0 | 265.6+15.0 
ARVD iPSC-CMs+3F | 4.340.6 | 44.0+4.3 | -32.0+3.0 | 348.0+18.9 
ARVD iPSC-CM+5F 39.24+3.7 | -33.9+2.5 | 317.6+13.9 


Table 1.1: The electrophysiological parameters measured in an iPSC human study 
of single cell electrophysiological ACM remodelling are summarised [46]. These cells 
present with downregulated subunits of the INa channel. The max upstroke velocity 
(Viner) and the peak action potential amplitude are governed by the sodium chan- 
nel. As we can see, PKP2 mutation alone (2nd row) and increasing levels of fatty 
infiltration (3rd and 4th rows) lead to a decrease in AP upstroke velocity and peak 
amplitude. The values of the minimum resting potential did not change significantly 
from the baseline values, nor did the corrected APD values at 90% repolarisation, 
indicating the predominantly a depolarisation disturbance. 


In ACM, the RYR2 receptor is shown to be the cause of leak of Ca?* in diastole, 
leading to Ca** depletion, arrhythmogenesis and heart failure or SCD [53]. What 
ends the calcium induced calcium release is an open debate. A few ideas are put 
forward, namely stochastic attrition, RYR2 adaptation, RYR2 inactivation and SR 
depletion [35]. Let’s look more closely at the mechanisms which lead to heart failure 
and result from a leaky RYR2 channel. 

Heart failure can be an advanced phenotype of ACM and can stem from inadequate 
cardiac pump function, which does not meet the metabolic requirements of the organs 
in terms of blood flow. Initially, failing cardiac function can be offset by increased 
beta adrenergic stimulation and catechomic signalling, whose mechanisms are outlined 
in Section 1.3.1. However, chronic adrenergic stimulation or prolonged stimulation 


during exercise are associated with defects in the functioning of the SR Ca?* release. 
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Figure 1.9: The differences between the Na* channel in controls (WT and Tg-WT) 
and over-expressed Dsg2 gene are shown in mice. The ion channel of the mutants 
shows a reduced peak negative potential (C) and reduced peak current (E) [50]. This 
gene mutation has been strongly associated with ACM development. An independent 
voltage clamp study of murine hearts confirmed these findings [51]. 


One such defective adaptation is a PKA hyperphosphorylation of the RYR2 channel 
[54] [55]. This makes the RYR2 channel exhibit a pathologically high level of open 
probability under resting conditions. Over time, this can lead to a reduction in the 
SR Ca?* content, with less available calcium and consequently weaker contractions. 

Above we have shown that a mechanism of spontaneous RYR2 opening in diastole 
can occur in particular physiological conditions. This results in a premature release 
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Figure 1.10: I-V relationship of Iy, in SCN5A mutant (red) and wild-type cells, 
derived from human embryonic kidney cells. Peak sodium is considerably reduced, 
INapeak = —27.7 + 2.7pA/pF (mutant) and Ina peak = —42.7 + 3.4pA/pF [52]. 


of SR Ca?", activating the forward mode of the sodium calcium exchanger. If the 
threshold to activate fast Nat is reached, an action potential will be triggered in 
this cell with the potential to lead to premature contraction of part or the whole 
ventricle [56] [57]. If the cell is in a region within the vulnerability window, this cell 
can become an ectopic or a re-entry circuit can be formed in its vicinity, leading to 
asynchronous activation of the ventricle. The same phenomenon was confirmed to 
lead to ventricular tachycardia in mice as a result of late RYR2 hyperphosphorylation 


during exercise [58]. 


1.4.1 Repolarisation abnormalities 


The Ix; channel has been found to be remodelled in ACM zebrafish, exhibiting a 
mutation in the plakoglobin JUP gene. The Ix; peak current is reduced by 30% as 
can be seen in Figure 1.11A, alongside a concommitant decrease in peak Jy, current 
(Figure 1.11B). 

cAMP levels are reportedly lower than controls in an ACM in vivo human study 
[59], suggesting an impaired beta-adrenergic response in ACM. Given that the beta- 
adrenergic cascade shortens the APD in control patients, it is expected that the APD 


will not shorten to as great an extent when ACM patients exercise. 
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Figure 1.11: Single cell ventricular zebrafish electrophysiology is recorded using patch 
clamping, in the context of a human mutation [2]. (a) Significant reduction in the 
Ix; zebrafish trace, shown on a current vs voltage plot. (b) In, peak current is 
concommitantly reduced in this mutation. 


Also known as PLN or PLB, phospholamban is a protein encoded by the PLN gene. 
It is a membrane protein, controlling the Ca?* pump in myocytes [60]. PLN is found 
to be a substrate for the cAMP-dependent kinase (PKA) in cardiac muscle. PLB 
is an inhibitor of the SERCA channel in the unphosphorylated state but inhibition 
ceases upon phosphorylation. Beta agonists limit the suppression of the SERCA 
pump so that Ca?* can be efficiently transported back into the SR during diastole. 
Studies show a mutation of the PLN gene is found in ACM patients, [61] suggesting 
ACM is also a disease of the calcium handling system under intense beta-adrenergic 
stimulation. 

A summary of key electrophysiological findings related to ACM and confirmed 
experimentally is included in Figure 1.12. 

Platelets 

Platelets, also known as thrombocytes, are a component of the blood which co- 
agulates blood vessel injuries. They do so by releasing a cocktail of substances, in 
the vicinity of the injured site. Some of these substances are believed to be arrhyth- 
mogenic and a comprehensive list of qualitative effects is summarised by Flores et al 
[62]. Recent dog [63] and human ex-vivo [64] studies have shown that platelet activ- 
ity is enhanced in ACM conditions. Due to their potential arrhythmogenic nature, 
platelet activity on electro-physiology (EP) needs further investigation, especially in 
the context of ACM. 

A list of the most physiologically important substances secreted by platelets is 


given here. 
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e ATP is released, partially inhibiting the Karp channel. [65]. ATP is also 


reported to increase transient and resting cytosolic Ca?*. 


e Platelet activating factor (PAF): In dog and pig studies, intracoronary injection 
reduces blood pressure, myocardial contractility and coronary blood flow [66], 
[67]. 


e Arachidonic acid: It is reported to disrupt calcium handling, increasing intra- 
cellular calcium and transient calcium influx across the sarcolemma. [68]. Pro- 
longed exposure to these acids leads to intracellular Ca?* overload in neonatal 


rat cardiomyocytes [69]. 


1.5 Effect of ACM gene mutations on structural heart 
muscle integrity 


Malicious genetic mutations have been discovered in desmosomes, a junctional net- 
work responsible for cell to cell adhesion [70] [71]. Genetics has suggested a new 
model for pathogenesis, hypothesising that impaired functioning of the cell junctions 
under stress leads to myocyte detachment and cell death, as well as inflammation. 
The body’s response is to provide a fibro-fatty replacement as repair. According to 
previous studies [6] [70] [72], at least 3 mechanisms contribute to the pro-arrhythmic 


substrate: 


e patches of myocarditis, or other local mechanisms, leading to cell death, fibrosis 


and cavities [73] 
e fibrous and adipose infiltrates within the myocardium 


e gap junction remodelling leading to structural disarrangement of the lattice of 


myocytes 


Several genes, responsible for the structural integrity of the ventricles, have been 
found to be related to ACM. The normal function of these genes regulates gap juction 
formation, cellular communication and cell coupling to the extracellular matrix. A 


summary of the sphere of influence of these genes is illustrated in Figure 1.13. 
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Figure 1.13: The cardiac desmosome and its proposed roles in (A) supporting the 
structural stability by ensuring cell-cell coupling (B) regulating the adipogenesis and 
apoptosis by controlling the transcription of the associated genes, maintaining smooth 
electrical conductivity through regulation of (C) gap junctions and (D) calcium home- 
ostasis. Figure adapted from previous genetics review on ACM [71]. 


1.5.0.1 Desmosomes 


Desmosomes couple adjacent myocytes together. This type of protein resists mechan- 
ical stress and can take 2 states, namely a strongly adhesive and hyper-adhesive state. 
The latter is characterised by a calcium independence, meaning that the mechanical 
properties of desmosomes are not influenced by the abundance of extra cellular cal- 
cium [74]. By contrast, the weaker state’s adhesion property is calcium dependent. 
Desmoglein-2 is a protein localised to the intercalated disk, responsible for mechan- 
ically and electrically coupling adjacent myocytes. It is encoded by the DSG2 gene, 
which is known to mutate in ACM. Studies on mice showed that mutations of this 
gene lead to lesions in muscles early on, followed by dying myocytes with calcification 
in the lesion areas. This led to widening of the inter cellular cleft, dissociation of 
intercalated discs around the lesions and a disturbed sarcomere structure. The dis- 
ease progresses into fibrosis and arrhythmias |75]. The main idea is that the mutated 
desmosomes are no longer able to support the mechanical stress of the contracting 
myocytes, leading to inflammation, induction of repair, fibrotic replacement, discon- 
nection of myocytes, permanent disarrangement of the myocardium and altered force 
distribution. The disease progression is outlined in Figure 1.14 and has a striking 


resemblance to human ACM progression [76]. Another Dsg2 mouse model of ACM 
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recapitulates the human phenotype of RV wall thinning, as a result of myocardial 
injury, extending transmurally and involving fibro-fatty repair and progressive disap- 
pearance of the myocardial muscle tissue |77]. 

DSP Gene 

Desmoplakin is a desmosomal protein encoded by the DSP gene. It enhances 
intracellular structural stability and its mutation is linked to cell disarray and fibro- 
fatty replacement in human variants of ACM [78]. A decrease in levels of desmoplakin 
and connexin-43 proteins at the intercalated disc was observed. A global inflamma- 
tion the ventricles, known as myocarditis, has been observed to be a symptom of 
ACM in 100% of patients in a cohort harboring a DSP mutation and may accelerate 
intracellular degeneration processes, myocyte death and fibro-fatty infiltration [79]. 

DSC2 Gene 

DSC2 encodes the Desmocollin-2 protein in cardiac myocytes. The protein is used 
in linking adjacent cells together in desmosome regions, much like the DSG2 gene. 
It is also shown to be involved in desmosome formation. Strong evidence suggests 
mutations are associated with ACM [80] [81], including mutations with a recessive 
inheritance [82]. Another mouse model found that a mutation in DSG2 leads to 
uncontrolled necrosis or muscle tissue disappearance [83] affecting both ventricles 
and involving concommitant fibro-fatty replacement. 

JUP Gene 

Plakoglobin is a protein encoded by the JUP gene. Plakoglobin is the cytoplasmic 
component of the desmosomes and adherens structures. Plakoglobin is essential for 
the development of intercalated discs and stability of cardiac muscle. More specif- 
ically, it is the only stabilising link between calcium dependent cadherins and the 
actin cytoskeleton. If the beta-catenin and plakoglobin proteins are suppressed in 
mice, this results in gap junction remodelling and arrhythmogenesis [84]. Addition- 
ally, dysfunction of plakoglobin results in a mutation of the heart progenitor cells, 
which could lead to adipocyte formation [85]. 

Connexins Connexins are a family of proteins forming gap junctions between 
cardiomyocytes. The connexin 43 (Cx43) protein is heavily expressed in the my- 
ocardial ventricular tissue and its abundance is associated with gap-junctional ion 
channel conductance levels, conduction velocity and structural gap junction stability 
[86]. Lack of this protein leads to significant dispersion of conduction, increased inci- 
dence of spontaneous arrhythmias, and loss of ventricular systolic function with only 


minor reductions in overall Cx43 expression in a murine model [87], particularly in 
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Figure 1.14: The hypothesis of disease progression in mouse as a result of DSG2 
mutation is outlined and figure is reported from [75]. Early phase is fully developed 
within 2 weeks of gene modification, acute phase phenotypes are visible by week 4 
and a much slower, chronic phase of the disease onsets from week 8 onwards. 


the context of the His Purkinje ventricular conduction system [86]. Cx43 is reportedly 
decreased in mutations associated with ACM [88]. 


1.5.0.2 Non-desmosomal mutations and structural manifestations 


AnkyrinG or AnkG 

AnkyrinG is a member of a family of adaptor proteins that mediate the attach- 
ment of the integral membrane proteins to the spectrin-actin based cytoskeleton [89]. 
AnkG is required for maintaining a polarised distribution of many membrane pro- 
teins including the Nat /K*ATPase Pump, the voltage gated Na* channel and the 
Na*/Ca?* exchanger. A physical association between AnkG and PkP2 has been 
observed recently [90]. A loss of AnkG expression leads to relocalisation of PKP2. 

PKP2 

PKP2 is a gene, whose mutation is known to be strongly associated with ACM [91] 
[92] and is one of the most frequently remodelled genes affecting cardiac structure, 
as shown in Figure 1.15. Mutations of this gene have repercussions on both electro- 
physiology and muscle structure, hence it is included herein. It encodes the protein 
Plakophilin-2, whose function is to link cadherins to intermediate filaments in the 
cytoskeleton, as can be seen in Figure 1.13. A loss of PKP2 expression leads to Iva 
remodelling [90], a decrease and redistribution of the connexin 43 protein as well as 
an overall decrease in coupling between adjacent myocytes. Furthermore, there is ex- 


perimental evidence in mice suggesting that plakophilin-2 binds to the ATP-sensitive 
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potassium channel, reducing its current density by 60% and introducing regional het- 
erogeneity in the channel [93]. Thus, the PKP2 gene and its associated protein can 
be thought of as a link between intracellular junction and cell excitability. 
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Figure 1.15: Gene mutation to macroscopic heart imbalances in ACM [94]. The gene 
mutation either influences the cardiac myocytes internally or externally. These effects 
combined lead to arrhythmogenesis. 


Silencing of PKP2 in cultured cardiomyocytes reduces the electrical conduction 
velocity [95]. It remains an open debate whether this reduction in velocity is due to 
reduced electrical coupling at gap junctions or decreased sodium current availability 
or due to both factors. 

The extra and intra-cellular mechanisms of cell binding and adhesion are explained 
in Figure 1.13. As the genes associated with structural integrity of the desmosomes 
mutate, they trigger a chain reaction, outlined in Figure 1.15. Either the gene mu- 
tation suppresses the intra-cellular process of cardiomyogenesis or it hits the extra- 
cellular network, comprising of ionic channels and gap junctions. The former leads 
to suppression of myogenesis, stimulating adipogenesis and the latter leads to cell 
death and decoupling/disarrangement of myocyte strands. A structural phenotype 


can include mild to severe, localised and global RV dilation [92]. The combination of 
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these effects results in arrhythmogenesis at heart level. 


% of 
Gene Locus Protein Disease 
PKP2 12p11 Plakophilin 2 25%—40% 
DSG2 18q12.1 Desmoglein 2 5%—10% 
DSP 6p24 Desmoplakin 2%—-12% 
DSC2 18q12.1 Desmocollin 2 o—T Vo 


Figure 1.16: An international consensus paper on ACM and other diseases outlines 4 
main genes are involved in the ACM mutation. PKP2 has the highest incidence [96]. 


PPARy 

Although not a gene, peroxisome proliferator-acitvated receptor gamma (PPAR y 
or PPARG) is a glitazone receptor, encoded by PARG gene in humans [97] involved in 
structural remodelling of the ACM heart. It is involved in the production of adipose 
tissue. Mice, which are genetically engineered to overexpress this gene, are found to 
develop dilated cardiomyopathy like symptoms [98], including heart dilation, reduced 
ejection fraction and increased levels of failure biomarkers. Human tissue sample as- 
says also found a link between PPRAGY, intracellular lipid overload, fatty production 
and contractile dysfunction [99]. Interestingly, this pathway is overexpressed in the 
RV yet is found to be functioning normally in the LV of patients studied, suggesting 
another mechanism may be at play in the LV fat generation process [100], which 
remains to be elucidated. 

SCN5A 

This gene has been shown to affect the structural muscle integrity in the heart, via 
its interaction with N-cadherin [52]. Depending on the protein interaction affected, 
genetic changes in NaV1.5 may thereby cause a predominantly structural phenotype 
(e.g. dilated cardiomyopathy), predominantly electrical phenotype (e.g. Brugada 
Syndrome, Long QT syndrome), or mixed electrical and structural phenotype (e.g. 
ACM). As such, the affected protein interaction, more than the exact gene mutation, 
may determine the phenotype. The biochemical mechanism governing the structural 
alterations due to SCN5A are yet unclear. This gene has originally been associated 
with control of the sodium channel [101]. Several mutations in this gene are observed 
to reduce peak Iya availability and prolong recovery from inactivation. This is one of 


the genes which impacts both muscle structure and single cell electrophysiology. 
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1.5.0.3 Non-genetic contributions to structural disease 


Intense physical exercise 

Competitive sport is observed to increase disease penetrance [102]. It promotes 
mechanical and oxidative stress. Areas of the ventricles show reduced reuptake of 
norepinephrine leading to chronic stimulation of adrenergic receptors, which in turn 
has been related to cardiac fibrosis. Fibrosis is known as the result of overactive 
fibroblasts, proteins which change the extracellular matrix properties of the tissue 
and result in a replacement of myocardial tissue into non-conducting scar [103]. It 
can result from inflammation, act as a reparative process when there is uncontrolled 
cell death, thus avoiding cardiac rupture or appear as a result of an overexpression of 
the TGF-6 gene. Fibrosis can take different shapes and affect heart muscle electrical 
and contractile function [104]. 

A major diagnostic criterion for ACM, RV dilation can be a result from intense 
physical exercise in mice with a PKP2 ACM mutation [105]. In this model of disease, 
connexin 43 delocalization at intercardiomyocyte gap junctions occurs in trained mice 


presenting with RV dilation. 


1.6 Mechanisms of conduction modulation in ACM 


There is an established interplay between conduction velocity and the balance of ions 
in cardiac muscle. There is an observed anisotropy of conduction in the longitudinal 
and transverse directions relative to myocyte strand orientations. Conduction veloc- 
ity is shown to be higher along the muscle fibers than across them [106]. Anisotropy of 
repolarisation has also been identified along the transmural direction of the ventricles 
and is of the order 20-25 ms [107]. This is physiological behaviour yet pathophysiolog- 
ical levels of the conduction velocity of the depolarisation wave in the ACM ventricles 
has been associated with low sodium channel conductance [108]. Conduction veloc- 
ity is modulated by the Jy, channel dynamics. It has been shown to be reduced in 
human and mice [50] tissue samples due to a lower density of sodium channels at 
the level of the intercalated disk. Mutations in the PKP2 and SCN5A genes associ- 
ated with ACM can lead to decreased peak sodium current, as discussed in Section 
1.4.1 of the thesis. An iPSC baseline model of ACM exhibited a 53% reduction in 
peak AP gradient upstroke, as observed in Table 1.1. This leads to a slowing of the 
electrical wave, as each cell will take longer to reach its peak potential and trigger 
the opening of the neighbouring cell’s voltage-gated Iy, channels. Patch clamp ex- 


periments in a mice PKP2 mutation of ACM [49], data from which are shown in 
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Figure 1.17, show a reduced peak Iyagz and Ig, currents. The qualitative effect of 
this mutation on AP morphology is an increased time to peak voltage and a less 
negative resting membrane potential. The combination of this electrophysiological 
remodelling on human ventricular tissue level conduction velocity is explored via a 
computational study by modulating Jy, and Ix; and studying the CV response in a 
fiber of human tissue, which can be found in Chapter 2. Moreover, in a human study 
on explanted ACM hearts [109], biological tissue analyses further showed disturbed 
distribution and significantly reduced expression of Connexin 43. A 50% reduction 
in this molecule has been recorded to severely impair conduction velocity in both the 
longitudinal and transverse fiber directions in mice [110] by tracking the electrical 
wave using an optical mapping technique, discussed in the imaging section of Chap- 
ter ??. These results are presented in Figure 1.18. The exact mechanism via which 
Cx43 influences conduction slowing remains unclear yet it is heavily expressed near 
gap junctional channels, suggesting it causes a disruptance to the passage of Nat ions 
between neighboring cells by affecting the conduction properties of the gap junctional 
ionic channels and delaying the electrical activation wave [111]. Connexin 40 and 43 
are heavily expressed in the fast conduction system of the ventricles, also known as 
the Purkinje network. This expression is known to increase conduction speed of the 
electrical wave in the Purkinje fibers [112] and a decrease in their abundance can lead 


to conduction slowing or conduction block in a recent mice study [113]. 
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Figure 1.17: A. 21 day old mice with a PKP2 mutation present with reduced peak late 
Inq current, when undergoing patch clamp I-V protocol. The data presented are the 
mean values for all the animals studied. B. The same experimental protocol revealed 
a decrease in peak current through the J, repolarising channel. C. The shape of 
the AP was altered as a result of these electrophysiological changes, inducing a less 
negative resting mebrane potential (RMP) during phase 4 of the AP in the PKP2 
mice model [49] . 
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Figure 1.18: A. Conduction velocity in the principal fiber direction is plotted against 
pacing frequency for control mice (white circles) and 50% induced Cx43 abundance 
mince (black). Notice that the percentage difference in CV increases for higher pacing 
rates. B. Conduction velocity in the tranverse muscle fiber direction is plotted against 
pacing frequency for control mice (white circles) and 50% induced Cx43 abundance 
mince (black). Notice that the percentage difference in CV increases for higher pacing 
rates also in this case [110]. 


In ACM, the number of gap junctions may be significantly decreased, while at the 
same time sodium current is down regulated. A study by Sato et al. [95] shows that 
the sodium current density can decrease by up to 50% in ACM human patients. It is 
not yet certain if there is a causal link between the two effects. 

As discussed earlier in this chapter, fibrosis may arise as the body’s healing re- 
sponse to myocardial degradation. It is another substrate which can alter the electrical 
properties of the myocardium. As a non-active conductor of electricity, the fibrotic 
scar’s electrical conductivity coefficient is lower than the surrounding myocardium 
[114]. The actual extent of conduction slowing may depend on the type of fibrotic 
infiltration, the extent as well as density of surviving myocardium within the scar 
region [115]. Fibrofatty replacement of the His-Purkinje system, affecting both its 
bundle branches, has been reported post-mortem in ACM [116]. This can lead to ma- 


jor conduction abnormalities of the ventricles as the Purkinje network is the highway 
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of conduction of the ventricles and contributes to uniform initiation of the ventricular 
activation sequence in the endocardium. Fat infiltration is another substrate likely 
to be a barrier to planar, smooth conduction. Fat is an insulator and is observed to 
produce conduction abnormalities in human patients imaged [117]. Cell disarray is 
also observed in human patient tissue samples and is associated with poor ventricular 
conduction [78]. A quantification of the contribution of these substrates to conduc- 
tion slowing is missing in the above mentioned studies and computer modelling of the 
ventricles can help to access the degree of conduction delay in regions where these 


substrates exist. 


1.7 Conclusion 


ACM is a rare form of potentially lethal inherited cardiomyopathy. Since the right- 
sided form of the disease was first defined in the 1980s, significant progress in under- 
standing the functioning of the disease has been made. Unfortunately, this has not 
led to successful umbrella treatments for the disease and patients have been treated 
using strategies developed for other forms of cardiac disease, such as anti-arrhythmic 
drugs, implanted defibrilators and lifestyle restrictions. Significant basic science ef- 
forts have contributed to undestanding the link between genes mutated and protein 
expression, with downstream effect on electrical and structural disease manifestation. 


A summary of key manifestations of ACM in the ventricles is given below. 


e Rupture of the intracellular myocyte coupling due to gap-junction remodelling 
leading to increase in RV chamber volume. 


e Fibrotic scars creation as a reparatory response to muscle degradation in the 


ventricular myocardium and in the fast Purkinje network of fibers. 


e Fatty infiltration as the body attempts to maintain the structural stability of 


the ventricles. 


e A decrease in peak sodium current depolarising the myocyte as well as decrease 
in peak action potential gradient dV/dt and a concommitant decrease in peak 


conductance of [x channel, leading to resting membrane depolarisation. 
e Contractility alterations as a result of Ca?* ion imbalances during the cell cycle. 


e Disarray of the fibers and cells within the muscle fibers. 
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e Myocardial muscle thinning due to uncontrolled cell death and consequent RV 


wall thinning. 


As can be seen, ACM affects the ventricles in several ways. The computational 
experiments performed for this thesis will aim to address how the electrical prop- 
agation during sinus rhythm in the ventricles changes as a result of ACM-induced 
remodelling. Shedding light on the dynamics involved will help understand disease 


function and increase our working understanding of this complex clinical entity. 
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Chapter 2 


Conduction Velocity Modulation due 
to Co-Expression of /y,—Ix, in ACM 


2.1 Abstract 


Conduction velocity (CV) heterogeneities are a long established pro-tachicardic sub- 
strate in cardiac disease, including ACM. One mechanism for generating dispersion 
of electrical conduction is conduction slowing in myocardial tissue. Both the Ij, and 
Ix ionic channels are known to have concommitant reductions of their conductances 
in given mutations of ACM in zebrafish [2] and human studies [49]. The Jy. channel 
is implicated in several ACM genetic mutations [118] [51]. The main role of [x1 in 
influencing CV may be via a co-modulation of the Iyq channel conductance yet this 
effect is yet to be quantified in human experiments. 

The rate dependent interactions of the Iv, — I ion channel complex on human 
ventricular CV are studied using an in-silico model. Two response domains of the 
myocardial tissue are found. There exists a slow frequency of pacing domain (1 Hz), 
where CV is largely governed by the Jy, channel alone. In the fast frequency domain 
(3.3 Hz), Ini plays a marginal role in determining ventricular CV. The simulation 
results suggest that, under arrhythmogenic conditions, the upper limit of the CV is 
governed by the Iva — [xn conductivity ratio. The simulations put an upper limit of 
45% conduction velocity modulation in ACM myocardial tissue, which is attributable 
to Ina — Ix, electrophysiological remodelling. Certain loop diuretics are known to 
induce severe hypokalemia in ACM patients [119] and this is suspected to lead to 
an increased count of arrhythmia. Another type of diuretic used in treating ACM, 
a potassium-sparing diuretic such as spironolactone, is known to prevent potassium 
excretion from the body and can lead to hyperkalemia if poorly managed [120] [121]. 
The role of the potassium extracellular concentration in promoting CV changes in 
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ACM ventricles is modelled and the findings discussed with respect to arrhythmia 


risk. 


2.2 Introduction 


Several cardiac diseases exhibit dispersion of electrical conduction velocity (CV), in- 
cluding atrial fibrillation [122], Brugada Syndrome [123], hypertrophic cardiomyopa- 
thy [124] and ACM [125]. Slow and heterogeneous CV, in combination with appro- 
priate anatomical substrates, translates into an increased risk of re-entry. Even in 
the absence of structural heart changes, such as in early-stage ACM patients with ion 
channel remodelling but otherwise healthy hearts [125], CV abnormalities have been 
shown to cause sudden cardiac death in young individuals [126]. Cardiac voltage- 
gated sodium channels (Jy) are transmembrane proteins responsible for regulating 
the flow of sodium ions. Mutations in the SCN5A gene, which reduce peak Jy, cur- 
rent in ACM [52], have been shown to lead to lethal arrhythmias [31]. The inward 
rectifier current (Ix) stabilises the resting membrane potential in phases 3 and 4 of 
the action potential (AP) cycle. [x1 has been experimentally shown to contribute to 
the creation of regional differences in cardiac excitability [127]. The role of Iva on 
governing CV in the context of ACM is well established experimentally. A study by 
Sato et al. [128] shows that the sodium current density can decrease by up to 50% 
in ACM human patients leading to alterations in CV. Similar results were obtained 
by optical voltage mapping in guinea pig ventricle, when altering the conductance of 
Tx, [129]. The co-regulation of Ivy, and Ig; is shown to occur in a plakoglobin ACM 
mutation, with data indicating a decrease in the conductance of both channels, as 
shown in Figure 2.1. Understanding the effect of co-modulating two important ionic 
currents, known to be responsible for CV modulation in isolation, is important in the 
context of ACM, as it can provide additional information on whether CV modulation 
is additive when both Iyq and Ix, are downregulated and help inform whether Iya 
and J, have similar quantitative effect on CV when co-modulated in human tissue 
simulations. 

Emerging experimental evidence suggests ionic channels are not regulated inde- 
pendently and are organised in macromolecular complexes [130]. The Iya-I«1 ionic 
channels form one such complex, which was shown to control cardiac excitability and 
arrhythmia by Milstein et al [1]. Their inter-species study concluded that an increase 
in functional expression of one channel reciprocally modulates the expression of the 


other, in order to maximise cardiac excitability. The same effect was observed when 
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the mechanism was one of reducing the maxiumum channel densities by infecting cells 
with an aderovirus. 

The extracellular potassium concentration [A *]° is another governing factor of CV 
in the myocardium. An elevated concentration of |A*]° , known as hyperkalemia, is 
experimentally shown to globally induce Jyq inactivation and [x1 opening [131], pro- 
moting longer refractoriness of myocytes. Equally, low levels of [A *]° are thought to 
lead to hyperpolarisation or an overall decrease in the resting cell potential. Both 
of these effects have repercussions on conduction. In the context of ACM, clinical 
reports have shown that ACM patients with structurally compromised hearts have 
more arrhythmic episodes in hypokalemic conditions [119].An investigation into the 
way hypokalemia affects CV can shed light on whether the arrhythmia generation 
mechanisms in ACM can be due to alterations in tissue CV. Hyperkalemia is exper- 
imentally known to induce arrhythmias, even without the presence of a pro-fibrotic 
structural substrate due to interaction with Jy, channel dynamics and APD short- 
ening [132]. Furthermore, hyperkalemia inducing kidney disease be found as a co- 
morbidity to ACM [133] and commonly prescribed diuretics such as spironolactone, 
used for the management of ACM symptoms [120] are known to increase [A*]°, can 
also lead to hyperkalemia. A quantification of CV modulation due to hyperkalemia 
can give order of magnitude estimates of the involvment of hyperkalemia in CV gov- 
ernance in the context of ACM patients. With this premise in mind, we set out to 
investigate the interplay between Iy, -[,, on ventricular CV and upstroke time in 
1D fibres. We focus on the rate dependence of the interaction by studying several 


pacing rates. 


2.3. Methods 


2.3.1 The Human Ventricular Cell and Tissue Model 


This in silico study was performed using a state of the art electrophysiological sin- 
gle cell model of the human ventricular cardiomyocyte, the Comprehensive in vitro 
Proarrhythmia Assay updated O’Hara Rudy dynamic (ORd) model [134]. It is an ex- 
tensively validated human cardiac electrophysiological model, incorporating the main 
ionic channels and sub-cellular processes involved in the myocyte AP cycle. It is used 
in the context of the bidomain formulation, in order to couple the cells in a tissue 
simulation. The bidomain formulation is a mathematical framework for describing 
the flow of charge within the myocardium. The assumption is that charge flows be- 


tween an extracellular and intracellular domain through ionic channels or due to a 
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Figure 2.1: Co-modulation of Iya-I«1 peak current densities in zebrafish cardiomy- 
ocytes when a plakoglobin gene ACM human mutation is introduced. A. Modulation 
of Ina peak current by (red trace) relative to a control Iyq black trace, in a voltage 
clamp experiment. B. Concommitant modulation of peak J, current in the same 
cell as in A (red trace), compared with a control zebrafish cell (black trace). 


potential gradient and can be described mathematically in Equation 2.1. A no flux 
boundary condition is imposed, meaning charge cannot leak out of the myocardium, 
as described by Equation 2.2. 


V.(a0;Vu;) + V(oeVve) = xCn( Se + Tion) (2.1) 


V.(a0V vo) = 0 (2.2) 


Where o; is the intracellular and o, the extracellular electrical conductivities, v; 
the intracellular potential, v, the extracellular potential, y is the the surface to volume 
ratio of the cell membrane, C’,, is the capacitance per unit area of the cell membrane 


and Jion the contribution from the cell ionic dynamics to the transmembrane current. 


2.3.2 In silico experimental setup 


Nodes are set up in 1D fibres, the simplest possible coupled system. The nodes are 
0.4mm in length and each node represents an ensamble of cardiac cells. The overall 
length of the fibre is 2e-m, comprising of 50 nodes. The model is paced into steady 
state initially before recordings take place, in order to ensure an equilibrium state of 
the model is reached. The 5 leftmost cells, situated at one end of the fibre, are excited 
simultaneously. Activation times (ms) are recorded by analysing the transmembrane 
potential, V, across each node. The activation time of the cell is defined as the 
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time point of maximum gradient of the transmembrane voltage. Activation times are 
recorded for a range of individual channel conductances, G [HF], via a parameter sweep 
for the 2 channels of this study, Iyq and I[x,;. The CV was subsequently calculated 
by the gradient of the least squares linear fit of the resulting plot of activation times 
of the cells. Simulations were conducted for a range of |0.1,1] scaling factors of the 
Inq and Ix, channel conductances, Gyg and Gx [4A], where 1 represents average 
current density and 0.1 is equivalent to 10% of the maximal current density. In order 
to measure if the Iy,_ - Ix, interaction is rate dependent, the fibres are paced using 
square waves of lms width, amplitude of -120000 pA/pF, at different periods or beat 
cycle lengths (BCLs). The 100 and 200 beats per minute cases are shown, where the 
first scenario mimics a resting heart rate in the upper limit of a control subject and the 
latter a fast ventricular tachychardia [135]. The healthy concentration of extracellular 
potassium ions in the model is [A *]° =5.4 mM. A scaling factor is applied to [K*]° 
in order to study the effect of abnormally high and low concentrations on ventricular 


CV. 


2.4 Results 


2.4.1 Low frequency pacing limit 


The CV at low pacing rates (> 600 BCL), was found to be almost entirely controlled 
by the Jy, conductance scaling factor. As one blocks a larger portion of the Inq 
channel, the CV also decreases, as illustrated in Figure 2.2. The CV varies by less than 
1% within the range of [;,; scaling factors studied, implying a weak dependence of the 
CV on Gx; in the low frequency pacing mode, equivalent to sinus rhythm. An ACM 
mutation is shown to reduce peak Iv, conductance by 84% and peak Ix; conductance 
by 71% [2] due to intercallated disk remodelling. By looking at the results in Figure 
2.2B, it can be concluded that an upper limit for conduction velocity modulation due 
to this effect is 45%, as the combined decrease in conductance of In,g-I¢ leads to a 
decrease of conduction velocity from 0.4m/s to 0.22 m/s. The above limit to ionic 
channel-caused CV modulation can be used to inform whole ventricular models of 
ACM, setting a lower limit for a possible ion channel-based CV slowing in myocardial 
tissue. The plots also suggest Ix, has a less dominant role in modulating CV in 
human tissue than murine models predict [129]. Using modulation factor values of 
Tna-L1 from an ACM zebrafish model with a plakoglobin mutation [2] to constrain 
the human ventricular analysis, suggests a decrease in tissue CV by up to 20% at 


slow pacing in human subjects, as per Figure 2.2(a). 


36 


CV [m/s] for BCL=600ms and control [K+]o 


eee 
0.2 0.4 0.6 0.8 1 


Scaling factor of GK1 


CV [m/s] for BCL=300ms and control [K+]o 


Ss 


o 
iv 
Scaling factor of GNa — 


° 
xu 


e 
x 


0.4 


° 
FS 


e 
a 
° 
a 


Cy 
z 
oO 
. 
3 
. 
6 
2 
ro] 
é 
a 
3 
a 
o 
a 


0.2 


0.4 0.6 0.8 1 
Scaling factor of GK1 


Figure 2.2: The CV along a 1D ventricular fiber [m/s] in case of (a) fast sinus rhythm 
pacing (600 ms BCL) and (b) tachycardic pacing (300 ms BCL) for different scaling 
factors of the original conductances, Gyg and Gx, . The CV is constant when varying 
the Ix, conductance and is governed by the sodium channel dynamics. The exper- 
imentally measured co-modulation of maximum current densities is plotted using a 
circle for a plakoglobin ACM mutation in a zebrafish model [2] and with a star for a 
rat model [1] not harbouring a specific ACM mutation. 


2.4.2 High frequency pacing limit 


The high frequency pacing range is not reachable in sinus rhythm and is presented 
here to simulate a tachicardic scenario. At fast pacing frequencies (BCL < 300 ms), 
the CV was found to be dependent on Jy - [x1 interactions. The CV exhibited a 
positive correlation with Gy, and a mild positive correlation with Gx, scaling factors, 
as illustrated in Figure 2.2(b), only at low values of Gx1<0.3. The Jj in isolation is 
found to not have a major effect on CV in human tissue simulations. These results 
are relevant in the context of ACM as some patients with a PKP-2 mutation [49] will 
still mostly have conduction velocity of the electrical wavefront depend on Jy, and 
only marginally on [x,. Yet, a co-modulation of Iy,q - In, may result in I, indirectly 
influencing CV by its modulating effect on Jy, . Using modulation factor values of 
Tna-Ix1 from an ACM zebrafish model [2] with a plakoglobin mutation to constrain 
the human ventricular analysis, suggests a decrease in tissue CV of up to 25% at fast 
pacing in human subjects, as per Figure 2.2 (b). Whilst experimental evidence for 
this process exists in murine non-ACM [136] and ACM specific zebrafish studies [2], it 
is yet to be quantified in human subjects. If this is the case, ACM mutations involving 
Ig, will have dangerous effects of reducing global CV and increasing heterogeneity of 
conduction as well as arrhythmogenesis in ACM patient hearts. 


2.4.3. High extracellular potassium reduces global CV 


ACM patients in advanced disease stages are prescribed spironolactone, a potassium 


sparing diuretic, which can cause hyperkalemia [137]. Furthermore, reports of ACM 
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patients developing kidney disease exist [133], which also has a consequence of in- 
creasing extracellular potassium levels. For this reason, the effects on tissue CV due 
to hyperkalemia are studied. Hyperkalemia features abnormally high extracellular 
concentrations of the potassium ion, denoted by [A*]°. It results from a high intake 
of [K*]° or kidney disease. It can be shown from the Nernst equation that an in- 
crease of extracellular A+ concentration raises the ion’s equilibrium potential. This 
is because the electrical potential opposes diffusion gradients across the cell mem- 
brane. Increasing [K*]° results in a smaller diffusion gradient and also a less negative 
equilibrium potential [138]. In turn, the resting membrane potential becomes less 
negative. As the recovery from inactivation of the [y, channel is voltage dependent, 
a more depolarised resting membrane potential inhibits a full sodium channel recov- 
ery. Furthermore, as In, ~ J/ [Kt]? in the ORd model, an increase in [K+]? implies 
a more abundant maximum J; repolarising current and hence a shortening of the 
APD. The same mechanism holds for the fast delayed repolarising potassium current 
Ixr, which is potentiated under hyperkalemia as its activity is dependent on [K*]?. 
The simulations confirm this mechanism, as shown in Figure 2.3. Hence, [vq peak 
availability is reduced and CV is slow in a hyperkalemic setting, as captured by Figure 
2.4. Qualitatively similar behavior was observed at a higher BCL of 600 ms. An ACM 
patient experiencing a hyperkalemic episode have further marginal reduction in CV, 
on top of any pre-existing conduction abnormalities due to substrate or ionic channel 
deficiencies. As the disease progresses with age, associated co-morbidities such as kid- 
ney disease, which has a hyperkalemic effect, have been reported in ACM [133]. This 
extra CV modulating factor will increase the likelihood of arrhythmogenesis also due 
to a shorteing of the action potential [124]. ACM patients are known to have a high 
ectopic burden|139] and an already slow and inhomogeneous ventricular propagation, 
hence the time window for ectopic-triggered reentrant arrhythmias will be even larger 


given a slower CV and a shorter action potential duration under hyperkalemia [140]. 


2.4.4 Low extracellular potassium illicits biphasic response of 
CV 


Hypokalemia, or an abnormally low concentration of [Kt]? in the blood is a well es- 
tablished pro-tachicardic factor. Recent animal studies show it can lead to prolongued 
ventricular repolarisation, slowed conduction and conduction block [141]. We have in- 
vestigated the effect of hypokalemia on cardiac electrophysiology and CV. Competing 


effects are observed and illustrated in Figure 2.4. 
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Figure 2.3: (a) A high extracellular [A *]° equivalent to 20% above control (6.48 mM) 
has the effect of hindering full recovery from inactivation of the Ivy, gate. This is the 
governing factor behind the slower CV observed. Notice that the resting membrane 
potential in (b) is higher due to the reduced diffusion gradient. 
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Figure 2.4: CV in a 1D fiber at BCL=300ms, at different concentrations of extra- 
cellular potassium. A scaling of 1 represents a control [K*]°=5.4mM. Notice that a 
hyperkalemic scenario rapidly decreases the CV for all combinations of the (Ga , 
Gx1) ionic complex studied. 
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The CV can be increased due to an increased Iyq availability due to an enhanced 
recovery from inactivation. A low [A*]° increases the diffusion gradient across the 
cell membrane and hence a more negative equilibrium potential is required to balance 
the flux of ions out of the cell, resulting in hyperpolarisation (Figure 2.5a). Contrary 
to the effect of hyperkalemia, hypokalemia fully opens the voltage-dependent Inq 
inactivation gates, once the cell has recovered (Figure 2.5b). Furthermore, the lower 
[A *]° availability directly reduces the function of the fast delayed potassium rectifier 
current Ix, (Figure 2.5c). This, in combination with a more prominent Iq spike 
and more negative resting membrane potential, increases the repolarisation time of 
the ventricular myocyte, leading to a prolonged APD, as experimentally confirmed 
in hypokalemia [141]. An ACM patient with already low CV in certain areas of the 
ventricle, owing to low Inq availability, connexin protein density and/or structural 
barriers to conduction is likely to experience a cardioprotective CV increase in the case 
of mild hypokalemia in the blood. Increasing conduction velocity is likely to reduce the 
risk of lethal arrhythmias [124], provided the increase in CV is globally homogeneous 
for the whole organ. This cardioprotective hypothesis can potentially explain thes 
reports of ACM patients recovering from persistent arrhythmias once their [A+]? 
values were raised back from severe hypokalemia (2.5mM) to approximately 4mM 
[119]. 

The system shows a biphasic response to low [K*]°. At very low [K+]? < 4.6mM, 
the delay in conduction due to lower resting membrane potential outweighs the speed 
up due to higher Jy, availability in our simulations, ultimately resulting in a slower 
CV for hypokalemic ventricular 1D fibers. This effect will be pro-tachicardic for the 
ventricles of an ACM patient [124], as it will increase heterogeneity of depolarisation 
times in the ventricles, which is shown to lead to increased arrhythmogenesis in 2D 
tissue simulations, when a structural scar substrate is already present [142]. The APD 
prolongation due to lower Ix, activity, mediated by a low [A*]° concentration is likely 
to be cardioprotective as it will reduce the number of possible re-entry combinations 
available to the wavefront as a larger part of the myocardium will be unexcitable at 
any one time. Confirmation of these hypotheses warrants a whole organ simulation 
study. 


2.5 Discussion 


Emerging electrophysiological consensus is that ion channels coexist in their macro- 


molecular multiprotein complexes [143] and alterations to genes and protein expres- 
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Figure 2.5: (a) The AP shape of a myocyte belonging to hypokalemic human ventric- 
ular single cell, parametrising [AK *]° = 3.78 mM or 30% less than the healthy control. 
Notice the higher peak amplitude of the hypokalemic simulation, due to an increased 
availability of Iva, as shown in (b), as well as prolonged repolarisation stage and 
lower resting membrane potential. The APD prolongation is due to the decreased 
Ix, current, whose output is negatively impacted throughout the cell cycle by the 
low [K*]° and shown in (c). 


sion leads to co-adaptation of the channels in the complex. It is important to un- 
derstand the mechanistic link between genetic mutation and an electrophysiological 


manifestation of disease in ACM, namely CV slowing, as it is one of the disease 


4] 


phenotypes present amongst many ACM genetic mutations in the ventricular human 
myocardium [118] [51]. In ACM, strong experimental evidence shows remodelling of 
the Iy,_ channel at single cell [144] and whole organ level [118]. Few studies have 
performed patch clamping and protein expression analysis on the J, channel yet 
data in zebrafish ACM models suggest combined modulation of Iya-[«K1, conductances 
[2]. Mean 1D fibre human ventricular CV was computed for a set of 3 parameters 
(INafast, Iki , [A *]° ). The Iva~ Ii ionic channel complex was confirmed to have 
an effect on global CV. This in silico study found that a positive correlation between 
Ix, conductance and CV velocity is present only at short cycle lengths (<300 ms) 
and is very marginal. At higher BCLs, the CV was found to be governed entirely 
by the Jy, channel conductance. A finding is the biphasic response of the CV un- 
der hypokalemic conditions is preserved in an updated version of the O’Hara Rudy 
model [134]. The original dependence of CV on [K+]? was shown computationally 
in a sheep Purkinje strand study [145], where the inverse parabolic CV profile with 
respect to [K*]° was first shown. The dependence of CV on Iya alone was also pre- 
viously studied computationally and the two were found to be intimately linked and 
that the L-type calcium current supports conduction only for values of Gy, below 
30% of control [146]. A speed up in CV is observed for mild hypokalemic episodes and 
a slow down of CV in acute episodes, indepedent of (Ina, Ix1) conductances. Based 
on experimental data in rat patch clamp experiments [1], the peak conductance co- 
modulation factor of Gy, when Gx, is reduced is 1.4, meaning a reduction of 50% 
in Gx, peak inward conductance will induce a 70% peak conductance reduction in 
Gna. Data on the reciprocal relationship was not found. Currently, there is no data 
in the literature to validate these findings in human cells, yet experimental work can 
pave the way for co-modulation being incorporated in computational human ventric- 
ular models. The work is contextualised with regards to ACM ventricular dynamics, 
namely that a reduction in Iv, conductance and/or Ix, via an indirect influence on 
Ina, & hyperkalemic or hypokalemic episode are all likely to increase susceptibility to 
life threatening arrhythmias in ACM patients given that they will reduce conduction 
velocity and increase heterogeneity of depolarisation times [147]. The only exception 
is in the case of mild hypokalemia (<10% change in baseline [A *]° ) , where the CV 


is expected to increase by a few percent. 
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Chapter 3 


MRI-informed modelling and 
simulation highlight the role of 
Purkinje and myocardial fibrosis in 
QRS abnormalities in 
Arrhythmogenic Cardiomyopathy 


3.1 Abstract 


The goal of this study is to investigate the electrophysiological and structural re- 
modelling explaining distinct QRS complex phenotypes in ACM. Human anatomical 
torso-biventricular models were constructed from cardiac MRI of ACM patients, with 
different degrees of QRS abnormalities. Computer simulations were conducted to dis- 
sect the contribution of ACM abnormalities such as 55% increase in RV end diastolic 
volume, 30% reduced wall thickness, conduction slowing in myocardium and in the 
Purkinje-endocardial layer. QRS complexes were reconstructed from the simulated si- 
nus rhythm activations and assessed based on their ability to reproduce ACM-specific 
QRS features, namely prolonged QRS > 100 ms, terminal activation duration (TAD) 
> 55 ms, low R-wave amplitudes < 2 mV, and loss of R-wave progression across 
precordial leads. 

Computer simulations demonstrate that, in ACM, mild QRS prolongation (+100 
ms), in absence of other QRS abnormalities, is explained by regional myocardial 
conduction slowing in the areas of MRI-identified fibrosis. However, MRI-based fi- 
brosis and RV chamber dilation alone proved insufficient to recapitulate more severe 
QRS abnormalities, such as severely prolonged mean QRS duration (+120 ms) across 
precordial leads, mean TAD across V1-V4 (+55 ms), and low QRS mean R-wave 
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amplitude (<2 mV) in precordial leads, as well as poor R-wave progression. A 70% 
decrease in Purkinje-endocardial conduction and 30% RV wall thinning were required 
to explain these more severe abnormalities. An analysis of additional substrates on 
the QRS complex shows the effects of global myocardial conduction slowing, RV 
myocardial conduction slowing, regional fibrosis-mediated conduction slowing and 
Purkinje-endocardial conduction system slowing all have distinct effects on the QRS 
complex features of an ACM patient. 


3.2 Introduction 


The diagnosis of ACM is challenging and requires complex scoring systems, combining 
multiple diagnostic modalities, namely the electrocardiogram (ECG) and magnetic 
resonance imaging (MRI) amongst others. Risk stratification is also complex due to 
the wide range of phenotypes affecting both ventricles and guidelines remain generic 
without addressing specific ACM phenotypes [148]. In the context of the evolving 
paradigm of precision medicine, identification of the pathophysiological mechanisms 
underlying different disease phenotypes in ACM is therefore critical to enable more 
patient-specific treatment recommendations and refine diagnosis. 

Diverse QRS abnormalities are clinically recognised in ACM, including prolonged 
QRS duration, QRS fragmentation, reduced QRS amplitudes, poor R-wave progres- 
sion, or prolonged S-wave terminal activation duration (TAD) in the absence of com- 
plete right bundle branch block [149]. Prolonged TAD has been included in the 2010 
Task Force criteria for diagnosis of ACM [3], and confirmed as one of the most dis- 
criminatory ECG features, used for differentiating ACM from other clinical entities 
in the latest 2019 expert consensus report [19]. Yet, the mechanisms underlying 
specific QRS abnormalities in ACM remain elusive. Several electrophysiological and 
structural components of the disease have been postulated as causal factors, including 
sodium channel dysfunction [144], right ventricle (RV) dilation [150], RV wall thinning 
[17], fibrofatty replacement of the myocardium [151] and/or of the fast conduction 
system [116]. Dissecting the contribution of each factor is challenging using clinical 
or animal studies, due to concomitant effects and current technical limitations. In 
particular, fibrofatty replacement, as identified by late gadolinium enhanced (LGE) 
MRI, is a diagnostic hallmark and collocalises with conduction slowing in ACM in 
ECGI studies [139]. 

Computer simulations using cardiac MRI-informed human torso-heart models al- 


low the quantitative assessment of the contribution of each electrophysiological and 
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structural factor on the ECG, as recently shown in hypertrophic cardiomyopathy 
[152]. Other in silico studies have shown the contribution of anatomical heart-torso 
variability on QRS morphology [153], the contribution of mechanical deformation on 
the ECG during systole [154], have assessed how acute ischemic events translate into 
deviations of the ST-segment and arrhythmic risk [155] and have investigated the 
contribution of heterogeneity in the activation sequence on the ECG in heart failure 
patients [156]. Whilst some computational studies investigating pump function in 
ACM exist [157] [158], no previous work has used in silico modelling and simulation 
to elucidate the origin of QRS abnormalities in this disease. 

The goal of this mechanistic study is to relate the key structural and electrophysi- 
ological ventricular remodelling to the QRS abnormalities identified in ACM patients 
using MRI-informed computer modelling and simulation. We quantify the implica- 
tions on the QRS complex of ACM-related remodelling, including RV dilation, wall 
thinning and fibrotic infiltration, as well as slow propagation in the myocardium and 
the Purkinje-endocardial layer, providing links between the underlying remodelling 
and the measured clinical QRS features. We hypothesise that, in ACM, myocardial 
conduction slowing in areas of LGE MRlI-identified fibrosis can explain QRS and 
TAD prolongation, without introducing changes in QRS morphology, whereas mor- 
phological QRS abnormalities such as low amplitude R-waves and abnormal R-wave 
progression, alongside severe QRS and TAD prolongation in the precordial ECG leads, 
are caused by abnormal Purkinje-endocardial activation and anatomical dilatation. 
Simulations using MRI-informed human torso-biventricular models highlight a promi- 
nent role of abnormal Purkinje-endocardial conduction in mediating morphological 
QRS interval abnormalities in ACM. 


3.3. Methods 


3.3.1 Clinical Data 


Twenty subjects, including 15 definite, 3 borderline and 2 possible ACM diagnosed 
patients following the 2010 ACM Task Force criteria [3] underwent a complete assess- 
ment at the Heart Hospital (London, UK), including high resolution, motion corrected 
1.5 Tesla cine MRI, LGE-MRI, ECGI, and genetic testing. Informed consent and local 
ethical approval were obtained for the study. LGE-/MRI based intensity thresholding 
was used to extract regions of fibrosis from the ventricles of the patients. A summary 
of the anatomical, structural and electrocardiographic evaluation of the cohort is pre- 
sented in Table 3.1. The data was published in full in a previous study [139]. 12-lead 
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ECGs were reconstructed from the body surface potentials, recorded in standard lead 
positions. The ECG features of a separate cohort of 42 definite ACM patients, ac- 
cording to the 1994 Task Force Diagnosis Criteria, were used to contextualise our 
ACM cohort [159]. 


3.3.2. MRI-based heart-torso models for ECG simulation 


MRI-based anatomical torso-biventricular models were reconstructed as illustrated 
in Figure 3.1, based on the MRI datasets [139], as described in more detail in the 
Appendix Section A.1.2 of this thesis. In brief, MRI images of the heart /torso were 
aligned to correct for breathing artefacts, and LV/RV endocardial, LV epicardial and 
torso surfaces were contoured. Initially, a uniform wall thickness of 3 mm was assigned 
to the RV, given the limitations of MRI to resolve the RV epicardial surface. The 
inability to resolve the RV wall thickness on all of the MRI stack is shown in the 
Appendix Figure A.2. Uncertainty in the location of the RV epicardial surface grows 
as the scanner moves from the RV base to the apex of the ventricle. The contours were 
used to fit 3D finite-element meshes of the ventricles and torsos of the subjects. The 
above geometry extraction pipeline was developped by Ernesto Zacur and adapted 
to ACM meshes by the candidate. The process involved testing and debugging the 
pipeline for new patient geometries, as well as working with the creator to devise 
algorithm and code improvements. 

Simulated activation sequences were computed on the constructed ventricular ge- 
ometries using the Eikonal equation [160]. This simplified model of electrical propaga- 
tion computes the time at which each microscopic muscle volume in the ventricles was 
excited, based on a shortest path algorithm. As in previous studies [152] [161] [162], 
the activation sequence was initiated using a phenomenological description of the ven- 
tricular trifascicular conduction system, by stimulating the inceptions (root nodes) of 
the free-running Purkinje system with the ventricular myocardium, then rapidly prop- 
agating on the endocardial surface. A conduction speed within physiological range 
(170 cm/s) was imposed on the endocardial layer, as calibrated in a previous study 
to obtain realistic QRS complexes [161]. LV locations of the root nodes were selected 
to correspond to the same anatomical locations producing a healthy reconstructed 
QRS sequence and validated against ex-vivo human activation sequence experiments 
[163]. A larger variability in the locations of the end points of the trifascicular con- 
duction system in the RV was observed by identifying the sites of earliest epicardial 
breakthrough on the ECGI epicardial surfaces. Hence, the root nodes in the thinner 
RV were adapted for each patient to obtain RV epicardial breakthroughs at the same 
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Figure 3.1: (left column) Computational pipeline for activation sequence and QRS 
complex simulations from MRI-informed torso/biventricular models was used to test 
hypotheses about the nature of remodelling in ACM patients’ ventricles. The ef- 
fect of anatomical (RV dilation and wall thinning), structural (myocardial LGE and 
Purkinje-endocardial layer fibrosis) and ionic (global and ventricle specific CV slow- 
ing) ACM remodelling on the QRS complex were studied. (right column) The clinical 
ECG data were post-processed and features of the signals were extracted, namely 
QRS, TAD and signal amplitudes. Patients were clustered in feature space. Features 
values and ECG morphology were compared between the clinical and simulation based 
ECGs on a patient basis in order to assess the validity of the substrate hypotheses in 
the simulations. 
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location as in the corresponding ECGI dataset. A more detailed description of the 
model parametrisation is provided in the Appendix Section A.2.2. The implementa- 
tion of the solution of the Eikonal equation was developed by Julia’ Camps, with the 
contribution to the method development of the candidate being the incorporation of 
fibrosis into the models and the distributed computing (parallelisation) setup of the 


code, which led to a 64 times speedup in computation time. 
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Figure 3.2: Modulation of conduction by fibrosis. A. Scaling of conduction velocity in 
LGE-regions for the three densities of fibrosis tested. No LGE-regions had a fibrotic 
intensity greater than 65% in our ACM patients, and hence conduction block was 
never modelled. B. Experimental modulation of conduction velocity in the presence 
of increasing fibroblast content [164]. Reproduced under CC-BY 4.0 Open Access 
license. 
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Once the activation sequence was simulated in the torso-biventricular anatomical 
model, the 12-lead ECG was computed at standard lead locations placed on the 
torso, as explained in Appendix Section A.1.3. In all figures, grid resolution for ECG 


visualisation is 40 ms/1 mV. 


3.3.3. MRI-informed modelling of ACM-related structural and 
electrophysiological abnormalities 


The effects of fibrosis, as identified by LGE MRI, on myocardial propagation were 
introduced by decreasing the conduction velocity (CV) in the patients in American 
Heart Association style regions of the ventricle affected (LGE+) [165]. For each 
subject, three different fibrotic density curves were used for their LGE-MRI informed 
models in order to account for uncertainty in conduction properties and given the 
large variability in fibrotic composition. If the region of fibrosis affliction extended up 
to the endocardial layer, CV in the fast propagation layer was equivalently reduced. 

Different models of fibrosis were tested in our study to account for the effect of 
varying degrees of fibrotic remodelling onto the 12-lead ECG, as illustrated by the dif- 
ferent curves of modulation of conduction velocity by fibrosis content shown in Figure 
3.2A. The shape of the tested curves was informed by established theoretical results 
on the square-root modulation of conduction velocity with respect to myocardial con- 
ductivity [166], as well as by experimental data on fibroblast-modulated conduction 
slowing [164], presented in Figure 3.2B. Hence, the uncertainty of conduction modu- 
lation via fibrotic scar tissue was explored by testing several strengths of fibrosis in 
our simulations and evaluating their impact on the QRS complex. To incorporate 
this information into the Eikonal model of conduction, the following approach was 
used. The Eikonal model is solved using a non-directed graph-based method [160], 
where the activation follows the shortest path through the graph. Each edge in the 
graph is originally assigned a weight based on the Euclidean distance between the two 
points it connects. An additional distance scaling factor was added between nodes in 
fibrotic regions, slowing the signal propagation between them. 

ACM is also known as a disease of the desmosomes [144], in which cell-cell decou- 
pling and other genetic causes can potentially lead to a slower global conduction in 
the myocardium. This was considered in the models by scaling the global tissue CV 
by factors in the range 25-95%. The results from tissue-level simulations in Chap- 
ter 2 suggest that up to 45% CV velocity modulation can be attributed to Ina-In1 
imbalances at cellular level. Further reductions can be attributed to structural gap 
junction loss of integrity, cell-cell uncoupling, changes in cell shape [4] and density 


a0) 


of channel forming connexin proteins [88]. Reductions of CV of more than 70% are 
thought to be unlikely in living patients yet are included to ellucidate the effect of 
CV modulation on QRS features more clearly. 

Furthermore, the effect of fibrofatty replacement affecting of the fast conduc- 
tion system [116] was investigated by considering up to 70% reduction in Purkinje- 
endocardial layer speed. This is supported by 3D electroanatomical mapping record- 
ings in ACM [159] showing that total endocardial activation times were significantly 
prolonged in ACM. 

Finally, ACM-related gene mutations lead to uncontrolled myocyte inflammation 
and subsequent apoptosis, where the myocardium is progressively replaced by a fi- 
brofatty substrate. This effect is especially prominent in the relatively thinner RV 
wall [167], where wall thinning can be very severe, with little surviving myocardium 
present at post-mortem examination. This effect was incorporated by reducing the 
thickness of the RV wall from 3 mm to 2 mm, in agreement with RV wall thickness 
reported for a cohort of 30 ACM patients [167]. The factors of ACM remodelling 
accounted for in this study are summarised in Figure 3.1b. 


3.3.4 ECG Features 


The simulated and clinical QRS complexes were compared focusing on ACM-related 
abnormalities, namely QRS duration, TAD, R-wave progression and R-S amplitude 
ratio, as highlighted in the clinical diagnosis guidelines [3]. The R-wave peak is defined 
as the maximum point on the curve when an R-wave is present. The Q-point on the 
precordial leads of the ECG is defined as the first point of negative gradient before 
the R-wave upstroke. The nadir of the S-wave is the minimum voltage of the QRS 
complex. TAD is defined as the mean across V1-V4 of the time interval from the nadir 
to the end of the S-wave, the point at which the voltage returns to the isoelectric value 
for the first time after the S-wave nadir [168]. QRS duration is the mean time elapsed 
between the Q wave start and S-wave end point across precordial leads. Clinically, 
a TAD of over 55 ms in the precordial leads and filtered QRS duration greater than 
114 ms are both minor ACM diagnosis criteria. The QRS features are labelled on 
the QRS complex in Figure 3.3A. A healthy R-wave progression across precordial 
leads is present when the R-wave amplitude grows steadily from V1 through the 
precordial leads, reaching a maximum before V6, then decreasing marginally. We 
propose the following criteria for the quantitative categorisation of the absence of 
R-wave progression (R+) in ACM: (i) low amplitude R-waves of less than 2 mV in all 
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precordial leads; (ii) R-wave to S-wave amplitude ratio less than 0.5 in V1-V4; and 


(iii) absence of inverse parabolic progression in R-wave amplitudes across V1-V6. 
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Figure 3.3: A. Example of QRS feature extraction. The QRS complex is defined 
as the time interval between the start of the Q wave (green circle) and the end of 
the ventricular activation, when the potential returns to the isoelectric value (black 
circle). TAD is measured from the nadir of the S-wave (turquoise circle) and the end 
of the QRS complex. B. QRS features in our ACM cohort do not cluster by genetic 
provenance (different colours). Patients with no genetic mutation were excluded from 


this plot. 


53 


Figure 3.3: C. Distribution of QRS features in different ACM cohorts. Open/filled 
blue symbols indicate LGE negative/positive (LGE-/LGE+) patients in our ACM 
cohort, respectively. Three patients did not exhibit evident S-waves and were ex- 
cluded from visualisation. Patients with definite (red circles) and borderline (black 
circles) loss of R-wave progression appear in the upper quadrant of protracted QRS 
width and TAD values. Results are presented for comparison against an independent 
ACM cohort (filled magenta symbols, no LGE data available). D. QRS complex in 
precordial leads of ACM LGE-+ patient with moderately prolonged QRS duration and 
TAD<55 ms (TAD-) and LGE transmural fibrosis in the LV anterior/anterolateral 
AHA segments. E. QRS complex in precordial leads of ACM LGE-+ patient with 
prolonged QRS duration of 120 ms and TAD>55 ms (TAD+), low amplitude signal 
(downward arrows), absence of R-wave progression and transmural LGE fibrosis in 
the RV anterior and anterolateral free wall in AHA plot. Grid resolution: 40 ms/1 
mV. Black lines highlight R-wave progression pattern. 


3.4 Results 


3.4.1 QRS-based phenotypes of ACM in LGE positive patients 


QRS duration versus TAD feature values for our ACM cohort are presented in Fig- 
ures 3.3B and 3.3C. Genetic mutations did not cluster in the feature space, and hence 
could not explain differences in feature values (Figure 3.3B). Regardless of LGE clas- 
sification, both QRS duration and TAD features progress from moderately prolonged 
to significantly protracted values, spanning the entire feature space, and highlighting 
an absence of correlation between LGE+/LGE- evaluation and QRS-derived features 
(Figure 3.3C, filled vs open blue symbols). The same distribution of QRS biomarkers 
was also observed in an independent ACM cohort [159], as shown in Figure 3.3C (ma- 
genta symbols). An absence of R-wave progression phenotype, spanning the entire 
range of the QRS duration and TAD feature space, was found in several of the ACM 
patients. Two very distinct QRS phenotypes, presented in Figures 3.3B-C, were taken 
from the feature space of LGE+ patients and studied in detail. 

On the one side of the spectrum is an ACM LGE-+ phenotype having no apparent 
QRS abnormalities. This phenotype presents a QRS duration of 100 ms across precor- 
dial leads, which is in the upper limit of the healthy range, and TAD<55 ms (TAD-). 
This was regardless of fibrotic replacement diagnosed by LGE-MRI, which was ob- 
served to coexist with delayed activation (high isochrones density) in fibrotic regions 
in ECGI examination [139]. Conversely, on the opposite side of the spectrum is an 
ACM LGE+ phenotype with multiple QRS abnormalities, involving prolonged QRS 
duration, TAD>55 ms (TAD+) as ACM Task Force 2010 criteria, and absence of R- 
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wave progression. The patient was LGE-+ in the RV anterior free wall, also revealing 
marked RV dilation and RV wall thinning in MRI examination. For completeness, the 
clinical presentations of both phenotypes will be referred to as LV_ LGE+/TAD-/R- 
and RV_LGE+/TAD+/R+, as outlined in Figures 3.3D-E. 


3.4.2 The role of right ventricular dilation on QRS features 


The role of RV dilation was examined by mapping an identical start of ventricular 
activation in positive and negative RV dilated patients, then allowing cardiac con- 
duction to proceed under healthy myocardial conditions, as detailed in the Appendix 
section A.2. 

The simulated activation sequences for the LV_LGE+/TAD-/R- and RV_LGE+/TAD+/R+ 
phenotypes are presented in Figures 3A and 3B, respectively demonstrating similar 
total activation times, except for delayed activation in the RV base of the dilated 
ventricle. Despite similarities in activation sequences, the inherently different ge- 
ometries resulted in substantially different QRS complexes (Figure 3.4C). The non- 
dilated ventricles gave rise to a normal QRS complex in all precordial leads, with 
regular R-wave progression (increasing from V1 to V5, and decreasing in V6), and 
QRS width and TAD in the normal range. On the contrary, activation on the dilated 
RV ventricles produced S-wave overshoots in V1-V3 and reduced R-wave amplitudes 
in V3-V4. This is explained by the dilation around the RV septal insertion point 
(Figure 3.4B, white arrows), and the fact that leads V3-V4 reflect septal activation. 
Hence, although RV dilation can produce significant QRS abnormalities in precordial 
leads, it was insufficient to explain the characteristic QRS features exhibited by the 
RV_LGE+/TAD+/R+ ACM phenotype. For reference, a simulated activation map 
and the associated QRS complexes of a control subject are included in Figure 3.4D 
and E respectively, illustrating a healthy activation sequence and R-wave progression, 


as no dilation or conduction abnormalities are present. 


3.4.3 LGE-MRI informed modelling explains clinical QRS fea- 
tures of the LV_ LGE+/TAD-/R- ACM phenotype 


We hypothesised that borderline QRS duration, in the upper limit of the normal 
range, is attributed to the late activating fibrotic region in the LV anterior wall, given a 
healthy ventricular anatomy in Figure 3.44. In order to model fibrofatty replacement, 
LGE-+ regions were imposed by introducing conduction slowing and testing different 


levels of conduction modulation, in order to account for uncertainty in the translation 
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Figure 3.4: Role of ventricular dilation on QRS complex. A. Control simulated acti- 
vation sequence for LV_ LGE+/TAD-/R- ACM phenotype, presenting no anatomical 
abnormalities. 
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Figure 3.4: B. Control simulated activation sequence for RV_LGE+/TAD+/R+ 
ACM phenotype, presenting marked RV dilation. Arrows indicate regions of preva- 
lent dilation at RV septal insertion. C. Reconstructed QRS complexes from results 
in panels A and B. The non-dilated geometry showed healthy QRS and R-wave pro- 
gression (left panel), while the dilated geometry (right panel) produced an abnormal 
S-wave overshoot and an absence of R-wave in leads V1-V3D. The activation sequence 
for a subject with control geometry. The root nodes are placed in standard locations, 
mimicking experimental findings by Durrer [163] for earliest epicardial breakthroughs. 
E. The simulated QRS complex precordial leads for the activation sequence shown in 
D , showing good R-wave progression across precordial leads with R-waves increas- 
ing from V1-V5 and decreasing in V6. . Isochrones separation: 10 ms. Black lines 
highlight the R-wave progression pattern. 


of LGE-MRI measurements into conduction properties. Simulated QRS complexes in 
precordial leads for low (40% slowing, light blue traces), intermediate (50% slowing, 
green traces) and high (60% slowing, orange traces) fibrotic densities are shown in 
Figure 44. LGE-MRI informed models exhibited increasingly prolonged TADs (and 
therefore QRS widths) for increasing values of conduction slowing, still preserving a 
normal R-wave progression up to a conduction slowing of 50%. If this threshold was 
exceeded (orange trace), QRS morphology deviated from the clinical phenotype. 
The corresponding epicardial views of ventricular activation for control, low, in- 
termediate and high fibrotic densities are presented in Figure 3.5C, exhibiting areas 
of increasingly dense isochrones at the LV base. Additionally, Figure 3.5D quanti- 
tatively summarises the derived QRS complex features for all the fibrotic densities 
investigated. Of note, the low density fibrotic region, located transmurally in the 
LV anterior/anterolateral free wall, correctly recapitulates a QRS duration of c.a. 
100 ms across precordial leads, whilst only mildly prolonging TAD in anterior leads 
V2-V4. Therefore, the LV_LGE+/TAD-/R- ACM phenotype could be explained by 
introducing LGE-MRI information into the model, and this directly translates into 
mild prolongation of QRS width and TAD in precordial ECG leads. Accessing the 
difference between the simulated and clinical QRS complexes as mean absolute error 
between pairs of signals was trialed. It was found to be a metric very dependent on 
the amplitudes of the signals and hence QRS feature analysis was the preferred way 


to access similarity between recorded and simulated QRS complexes. 
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Figure 3.5: A. Effects of fibrofatty replacement, present transmurally in the LV an- 
terior /anterolateral segments, in the LV_ LGE+/TAD-/R- ACM phenotype. 
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Figure 3.5: A. QRS complexes in precordial leads for control versus low (light blue 
trace), intermediate (green trace) and high (orange trace) density fibrotic patches, 
located in the LV anterior and anterolateral free wall. Only the most dense fibrotic 
patch leads to loss of R-wave progression. Black lines highlight the R-wave progression 
pattern. The mean absolute error is with respect to the clinical QRS. B. Clinical 
QRS complex of the LV_ LGE+/TAD-/R- patient. C. Activation maps for the same 
conditions as in A. D. Quantification of QRS duration and TAD for different fibrotic 
densities in LGE-+ regions. The difference between the simulated (A) and clinical (B) 
QRS complexes was reported as mean absolute error between pairs of signals. 


3.5 LGE-MRI informed modelling alone is insuffi- 
cient to reproduce the QRS morphology of the 
RV_LGE+/TAD+/R+ ACM phenotype 


Next, we assessed whether LGE-MRI information is also sufficient to explain the 
RV_LGE+/TAD+/R+ ACM phenotype, in terms of reproducing its characteristic, 
extremely prolonged QRS width in all precordial leads, protracted TAD in V1-V4 
(including lateral leads V5-V6), low amplitude R-waves, and abnormal R-wave pro- 
gression. 

We translated the LGE+ information, indicating alarge transmural fibrotic patch 
in the anterior/anterolateral RV, into the personalised anatomical model, following 
the same methodology as for the LV_ LGE+/TAD-/R- ACM phenotype. In addition, 
we personalised the Purkinje inceptions on the RV endocardium in order to replicate 
the epicardial breakthroughs observed on ECGI. Introducing LGE-MRI information 
into the model delays activation in the RV anterior and lateral free walls, as high- 
lighted by the activation sequences in Figure 3.6C, eliminating the S-waves in leads 
V1-V4 (Figure 3.6A). This significantly deviates from the clinical phenotype (Figure 
3.6B), where wide and deep S-waves in leads V1-V4 are observed. Hence, LGE-MRI 
information alone is not sufficient to explain the QRS abnormalities of this phenotype. 

The capacity of dysfunctional sodium channels to account for the QRS abnormal- 
ities observed in this ACM phenotype was also tested by reducing global conduction 
speed in the LV and/or RV, as well as in isolation or concomitantly to fibrofatty 
replacement [144] [101]. No combination of global conduction slowing and/or fibro- 
sis could account for the QRS features of the RV_LGE+/TAD+/R+ ACM patient 
phenotype, as summarised in Figure . 
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Figure 3.6: Effects of fibrofatty replacement in the RV_LGE+/TAD+/R+ ACM 
phenotype. A. Simulated QRS complexes for control versus low, intermediate and 
high density fibrotic patches, located in the anterior and anterolateral RV. The black 
lines indicate the R-wave progression pattern. The mean absolute error is given with 
respect to the clinical QRS. 
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Figure 3.6: B. The clinical QRS complex for the RV_LGE+/TAD+/R+ phenotype. 
C. Activation time maps for the same conditions as in A. In the absence of additional 
factors, LGE+ information alone produces QRS complexes which deviate from the 
clinical phenotype in leads V1-V4. The difference between the simulated (A) and 
clinical (B) QRS complexes was reported as mean absolute error between pairs of 
signals. 


3.5.1 Fuibrofatty replacement of the fast conduction system ex- 
plains prolonged TAD and absence of R-wave progres- 
sion of the RV _ LGE+/TAD+/R+ ACM phenotype 


Fibrofatty replacement of the His-Purkinje system, affecting both its bundle branches, 
has also been frequently reported post-mortem in ACM [116]. However, it remains 
overlooked as a causal mechanism of the disease, possibly due to current limitations 
in LGE-MRI resolution and lack of noninvasive assessment. 

In order to ascertain the repercussions of fibrofatty replacement in the fast conduc- 
tion system on QRS features, the Purkinje fast endocardial layer conduction speed 
was reduced in both ventricles, mimicking cellular uncoupling and fibrosis in the 
Purkinje-endocardial layer, while keeping the bulk tissue conduction speed unaltered. 
The resulting QRS complexes and activation sequences for a range of Purkinje con- 
ductions speeds are shown in Figures 3.7A and 3.7C, respectively. A reduction in R- 
wave amplitude was identified as a result of the decrease in the Purkinje-endocardial 
layer speed, together with increases in QRS duration and TAD. The clinical phe- 
notype of R-wave progression loss (Figure 3.7B) was reproduced by a 70% slowing 
in the Purkinje-endocardial layer (Figure 3.7A, green trace). A functional Purkinje- 
endocardial layer results in an almost instantaneous activation of the endocardium, 
followed by transmural propagation, and finally activation of the epicardium (Figure 
3.7C, control). Slowing the Purkinje-endocardial layer results in a slowed activation 
of both the endocardium and epicardium (Figure 3.7C, slow Purkinje CV). In partic- 
ular, a marked slow conduction in the Purkinje-endocardial system reproduced most 
of the QRS abnormalities exhibited by the RV_LGE+/TAD+/R+ ACM phenotype, 
namely QRS duration of over 120 ms, TAD>55 ms across leads V1-V4, and absence 
of R-waves in leads V1-V4, yet was unable to match the reduced R-wave to S-wave 


amplitude ratio observed in the precordial leads QRS complexes of the clinical ECG. 
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A Simulated QRS Complexes 
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Figure 3.7: Effects of fibrofatty replacement of the His-Purkinje system in the 
RV_LGE+/TAD+/R+ ACM phenotype. 
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Figure 3.7: A Simulated QRS complexes in precordial leads for varying Purkinje 
conduction speeds (control: 0% reduction; intermediate: 35% reduction; slow: 70% 
reduction) in an otherwise healthy myocardium. The orange trace shows how a reduc- 
tion in RV wall thickness translates into sharper R-waves, deeper S-waves in V1-V4 
and smaller TAD with respect to the slow Purkinje QRS. The black lines show R- 
wave progression. B. The clinical QRS complex for the RV_LGE+/TAD+/R+ phe- 
notype. C. Activation time maps for the same conditions as in A. D. QRS features 
of the models in A. Marked slow conduction in the Purkinje system recapitulates the 
clinical phenotype, including low R-wave amplitudes, absence of R-wave progression, 
prolonged QRS width and TAD in leads V1-V4, and a prolonged QRS duration in 
V5-V6. The difference between the simulated (A) and clinical (B) QRS complexes 
was reported as mean absolute error between pairs of signals. 


3.5.2 RV wall myocardial thinning produces deeper S-waves 
on the anterior /anterolateral leads 


The only remaining QRS feature of the RV_ LGE+/TAD+/R+ ACM phenotype un- 
explained by a combination of anatomical and fibrotic remodelling is the observed 
reduced ratio of R to S-wave amplitudes in anterior/anterolateral leads. Quantita- 
tively, the mean R-wave to S-wave amplitude ratio in leads V1-V4 was 0.42 in the 
clinical ECG, whilst for the slow conduction Purkinje scenario (Figure 3.7A, green 
traces) it remained 0.58, suggesting less deep S-wave deflections than in the clinical 
recordings. 

The RV wall thinning observed in the clinical evaluation of the patient was incor- 
porated at this point in the ventricular reconstruction. A reduction of the RV wall 
thickness from 3 to 2 mm resulted in significantly deeper simulated S-waves and the 
subsequent decrease in R-wave to S-wave amplitude ratio to 0.37, otherwise preserv- 
ing the rest of QRS complex features, as shown in Figure 6A (orange trace). This is 
explained as follows. The transmural activation of the RV wall translates into positive 
deflections on the ECG as the depolarisation front travels towards the recording elec- 
trodes from the endocardial to epicardial layers, producing less prominent R-waves 
in the thinner myocardium [168]. Once the activation reaches the epicardial surface, 
it translates into negative deflections as the depolarisation front travels away from 
the recording electrode. The thicker the RV wall, the longer the transmural depo- 
larisation, effectively contributing to reduce the overall S-wave amplitude because 
transmural wave propagation will be occurring at the same time as epicardial prop- 
agation in different parts of the ventricle. This secondary opposing effect is thus 
reduced when the RV is thinner, leading to a deeper S-wave morphology. A thin RV 
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resulted in a shorter QRS complex when the Purkinje conduction was slowed (Figure 
3.7D, orange column quantifying QRS duration). The QRS features quantification 
showed that both the 70% reduced Purkinje-endocardial slowing alone (green bar) 
and in combination with RV wall thinning (orange bar) lie close to the clinical QRS 
features (red bar) in Figure 3.7D. The QRS complex duration produced with the slow 
Purkinje-endocardial layer and thin RV (Figure 3.7 D, orange trace) is reduced in V1, 
V2 and V4 with respect to when only a slow Purkinje-endocardial layer is considered. 
This is because some of the latest activating regions appear in the base of the right 
ventricle, hence reducing the thickness of the RV slightly shortens the total activation 


time. 


3.6 Effect of CV slowing substrates on QRS features 


In order to assess the role of each type of remodelling considered on the clinical 
QRS biomarkers, we compute the Spearman’s partial correlation coefficient between 
model inputs and outputs. By performing a grid search over the degree of Purkinje- 
endocardial coupling conduction velocity, myocardial LGE scar velocity, whole my- 
ocardial tissue conduction velocity and RV tissue conduction velocity and quantifying 
the precordial QRS complex reconstructed from the simulations, it is possible to infer 
what is the isolated effect of each of these substrates on the QRS complex and TAD, 
as well as on the R. to S wave amplitudes ratio. The RV was considered as an isolated 
target for remodelling, as we hypothesise that the anatomical remodelling in the thin 
RV could translate into cellular uncoupling and reduced conduction velocity due to 
gap junction dysfunction. Differences in sodium channel densities between the two 
ventricles have also been highlighted in recent canine studies [169]. 

The correlations are summarised in Figure 3.8A. When analysing a control simula- 
tion without CV slowing substrates for patient RV_LGE+/TAD+/R-+, a mean QRS 
complex duration of 58.9ms across precordial leads is simulated with mean TAD of 
less than 29.3ms and a mean R-S wave ratio in leads V1-V4 of 0.91. The clinical mean 
QRS complex duration is 121.3ms across precordial leads, mean TAD is 58.0ms and 
R-S ratio is 0.43 across V1-V4. A type or combination of remodelling is sought to con- 
commitantly increase QRS complex and TAD, whilst decreasing R to S-wave ampli- 
tudes ratio, in order to reproduce the clinical features of the RV_LGE+/TAD+/R+ 
patient. The simulation results show that reducing myocardial bulk tissue conduc- 
tion velocity (CV) increases both QRS duration and TAD (Figure 3.8B(left panel)), 
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yet also increases R-to-S-wave ratio (Figure 3.8C(right panel)). This therefore de- 
viates from the investigated RV_LGE+/TAD+/R+ phenotype (characterised by a 
low R-to-S-wave amplitude ratio). A reduction in the Purkinje-endocardial layer CV 
prolongs both the QRS and S-wave TAD duration of the precordial leads, whilst de- 
creasing the R-to-S-wave amplitude ratio. Introducing either LGE-informed or global 
RV conduction slowing increases QRS duration, yet decreases TAD and increases R- 
to-S-wave amplitude ratio, deviating from the clinical RV_LGE+/TAD+/R+ QRS 
phenotype. 

In summary, all the types of CV modulating remodelling tested increase QRS du- 
ration when a reduction in CV was applied to this specific patient, yet only a reduction 
of the Purkinje-endocardial CV concomitantly increases TAD and decreases R-to-S 
wave amplitude ratio. Hence, only the reduction of Purkinje-endocardial layer CV 
at control levels of bulk myocardial CV could reproduce the RV_LGE+/TAD+/R+ 
QRS phenotype. 


3.7 Discussion 


ACM is a complex disease involving desmosomal and gap junction uncoupling, sodium 
channel downregulation, followed by fibrofatty replacement, leading to conduction 
delays in the myocardium. The mechanistic basis of the QRS features in this disease 
remain unclear, especially the factors determining different QRS phenotypes in ACM. 
This translational modelling study, combining anatomical and structural data from 
LGE-MRI, explores the mechanistic link between tissue pathophysiology and clinical 
QRS phenotypes. The main findings of this study are: 


e The LV LGE+ ACM phenotype can be fully explained by introducing my- 
ocardial fibrosis as informed by LGE-MRI into personalised models of ACM 
ventricular anatomy. This directly translates into the upper limit of normal 
QRS duration and TAD in precordial ECG leads. 


e The RV LGE+ ACM phenotype, on the contrary, is not explained by myocardial 
LGE-MRI information alone. This more complex phenotype, characterised on 
the ECG by severely prolonged QRS width, prolonged TAD, and absence of 
R-wave progression, further requires a model presenting with fibrotic Purkinje- 
endocardial conduction system and RV wall thinning to recapitulate its clinical 


signatures. 
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Figure 3.8: Correlation between remodelling and QRS biomarkers. A. Pearson’s 
correlation coefficients between different CV scaling factors and QRS biomarkers. 
CV was scaled by the factors [0.4, 0.6, 0.8, 1] in each of the types of remodelling, 
giving rise to a total of 256 bi-ventricular models considered. B. Change in mean QRS 
duration and mean TAD across precordial leads (left panel), and R-wave to S-wave 
amplitude ratio (right panel) across leads V1-V4. Results are shown for slowing of 
the bulk tissue CV in the RV_LGE+/TAD+/R+ ACM ventricular geometry, while 
keeping the rest of CV scaling factors constant to 1. The left panel shows the increase 
in QRS duration and TAD as bulk myocardial tissue CV is decreased. The right 
panel shows that R-to-S-wave amplitude ratio increases when global myocardial CV 
velocity is decreased, deviating from the RV_LGE+/TAD+/R- clinical phenotype, 
which has a lower R-to-S-wave amplitude ratio than the control model. 


Conduction abnormalities in ACM are well reflected in the QRS complex, includ- 
ing alterations in its duration, morphology or amplitude [170]. In particular, QRS 
duration is markedly prolonged in ACM especially in V1-V3, in up to 64% of patients 
[171]. Within the QRS complex, a protracted TAD has received special attention 
for diagnosis, as it is present in 47-71% of ACM patients [172] and was confirmed 
as highly discriminating from other clinical entities [19]. Indeed, both the 2010 Task 
Force ACM diagnostic criteria [3] and the latest 2019 ACM expert consensus state- 
ment [19] consider a TAD>55 ms as an established ECG diagnostic criteria of ACM. 
In our results, presence of fibrosis directly translated into moderately protracted QRS 
duration and TAD in precordial ECG leads. 

Poor R-wave progression in precordial leads, in the absence of bundle branch block, 
is found in 22% of ACM patients [173], and is highlighted as one of the hallmarks of 
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the RV_LGE+/TAD+/R+ phenotype. In this study, we introduce quantitative cri- 
teria for diagnosing the absence of R wave progression (R+ phenotype) in ACM, and 
observe that span the phenotype spectrum of QRS width and TAD values. Our results 
highlight that a slow Purkinje-endocardial conduction in ACM leads to the loss of 
precordial R-wave progression. The remodelling of the fast Purkinje-endocardial layer 
was indeed the only mechanistic factor able to replicate the loss of R-wave progression, 
low R-wave amplitudes, and extremely prolonged QRS duration and TAD, character- 
ising the RV_LGE+/TAD+/R-+ phenotype. Importantly, the Purkinje conduction 
slowing was found to affect both ventricles equally in the RV_LGE+/TAD+/R+ 
phenotype. This remodelling is in agreement with previous clinical studies on en- 
docardial activation, suggesting severely prolonged total RV endocardial activation 
in ACM patients, as measured by 3D electroanatomical mapping [159]. In a sepa- 
rate study, endocardial activation was delayed by 70% with respect to the mean of 
controls in an ACM cohort without significant structural disease [125], suggesting 
that Purkinje-endocardial remodelling may precede myocardial fibrosis. Moreover, a 
dysfunction of the His-Purkinje system was reported in 68% of the hearts (54% with 
fibrous, 8% with adipose and 5.6% fibrofatty infiltration) in post-mortem studies of 
200 ACM-related sudden cardiac deaths [116], which could also explain the endocar- 
dial conduction slowing in the electroanatomical mapping studies. In fact, the authors 
of that study note that Purkinje dysfunction in ACM may have been overlooked and 
that two thirds of deaths involve Purkinje remodelling, while the main focus so far 
has been on the RV wall. Different genetic mutations in ACM have been also shown 
to reduce the expression of connexins 40 and 43 [174] [175], both highly expressed in 
the His-Purkinje conduction system. This may provide a molecular explanation of 
the reduced Purkinje-endocardial conduction identified in our study as causal factor 
for this phenotype. 

Finally, RV wall thinning reproduced the low R-wave to S-wave amplitude ratio 
observed in the RV_LGE+/TAD+/R+ phenotype at the expense of reducing QRS 
duration. This may be a confounding factor at diagnosis, as patients with increased 
disease severity, including RV wall thinning, may present with only mild QRS pro- 
longation. 
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Chapter 4 


Conclusions 


This study set out to investigate the mechanisms of disease underlying Arrhyth- 
mogenic Cardiomyopathy with the idea of furthering mechanistic understanding of 
disease function using computational and imaging techniques. A vast literature re- 
view revealed that few computational studies had investigated ACM and that even 
going beyond the computational work and into clinical or basic science approaches, 
ACM remains an understudied entity. Diagnosis is challenging in this clinical entity, 
with a complex scoring system being employed to diagnose patients. The 12-lead 
ECG is heavily relied upon for diagnosing ACM and several other diseases such as 
Brugada Syndrome, Hypertrophic Cardiomyopathy and other Cardiomyopathies have 
overlapping clinical presentation [176]. Risk stratification is even more difficult, as 
the disease involves several electrophysiological, structural and geometrical substrates 
with overlapping features on the QRS complex. The study set out to quantify the 
contribution of ventricular remodelling to organ level conduction abnormalities and 
their consequent manifestation on the 12-lead ECG patterns. Ventricular types of 
remodelling considered in the simulation study include combined effects of channel 
conductance alterations Iy, — Ix, on global ventricular conduction velocity, geomet- 
rical effects such as RV wall thinning due to uncontrolled apoptosis and RV enlarg- 
ment on activation sequence dynamics and fibrosis-mediated conduction slowing of 
the Purkinje-endocardial fast conduction system and of the bulk myocardial tissue. 
Preliminary simulations in 1D strands setup found that a co-modulation of Ia — Iki 
in ACM could reduce conduction velocity in tissue by up to 45%. Additional decreases 
in conduction velocity would therefore be the result of additional coexisting types of 
remodelling. Patients were assessed in terms of their QRS complex features, namely 
magnitudes of QRS, TAD and state of R-wave progression across precordial leads. 
Two distinct ACM patient phenotypes were identified in the medical cohort and their 
QRS complex features were explained in terms of combinations of different types of 
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ventricular remodelling. This was achieved by running simulations of sinus rhythm 
ventricular activation sequence on a personalised ventricular anatomy extracted from 
MRI images of the patient’s ventricles. 

The study found that a patient presenting with a long yet physiological QRS 
complex duration (<100ms) and normal values of TAD (<55ms), as well as with an 
LGE-MRI localised fibrotic scar in the anterior left ventricular wall was recapitulated 
by in-silico simulations when introducing 40% conduction slowing in the LV area 
presenting with fibrotic remodelling and personalising the contact points of the fast 
running Purkinje tree on the RV endocardial surface. The fibrotic-scar mediated con- 
duction velocity slowing resulted in an increase of the precordial leads QRS complex 
duration by 25ms and a prolongation in TAD in leads V2-V4 by a similar amount. 
A more advanced QRS phenotype involving loss of R-wave progression, extremely 
protracted QRS duration (>120ms) across precordial leads, prolonged TAD (>55ms) 
across V1-V4 and low R to S-wave amplitudes ratio could not be reproduced in the 
simulations by introducing only an LGE-informed fibrotic scar, present in the RV 
anterior wall and RV chamber dilation in this patient. Fibrotic infiltration-mediated 
conduction slowing of 70% in the fast Purkinje-endocardial conduction layer, along- 
side 30% RV wall thinning and RV chamber dilation were the substrates required 
to produce the aforementioned QRS abnormalities in the simulated QRS complexes. 
Furthermore, the effects on the QRS features as a result of introducing several dif- 
ferent types of remodelling were studied using the advanced patient geometry. It 
was found that reducing global myocardial conduction velocity increased QRS and 
TAD values yet also increased R to S-wave amplitudes ratio in precordial ECG leads. 
Decreasing conduction velocity only in the RV led to increases in QRS duration and 
R to S-wave amplitudes ratio but decreased TAD as it resulted in progressive disap- 
pearance of the terminal S-wave. Decreasing conduction in the fibrotic patch in the 
RV anterior wall was found to increase QRS complex yet also led to disappearance of 
the S-wave and hence a reduction in TAD and an increase in R to S-wave amplitudes 
ratio. A global conduction velocity reduction mimicking cell-cell gap juctional remod- 
elling and low Ing — Ix availability was found to increase QRS complex, TAD and 
R-S wave amplitudes ratio. Merely introducing RV chamber dilation alone resulted 
in loss of R-wave progression across precordial leads and a complete loss of R-wave 
in leads V1-V3. 

In summary, graph-based personalised ventricular simulations of electrical activa- 
tion were used to explain two distinct ACM precordial QRS complex phenotypes. The 
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role of individual and combined tissue level ventricular substrates on the QRS com- 
plex features was explored. The computational pipeline presented herein can be used 
to refine diagnosis non-invasively and quantify the extent of ventricular remodelling 
present in the ACM ventricles of patients. The work can be embedded into clinical 
pipelines to stratify risk and guide therapy in ACM and adds confidence in the added 


value provided by the use of computer simulation techniques in bio-medicine. 
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Appendix A 


Appendix A 


A.1 MRI Processing 


A.1.1 Construction of MRI-based biventricular meshes 


Recently, we presented a computational pipeline for the generation of personalised 
torso-biventricular meshes from clinical MRI datasets [153] [177]. The first step is 
to create a surface ventricular mesh from sparse, non-parallel, cross-sectional, non- 
coincidental, endocardial and epicardial contours, as published in our previous work 
[153]. The approach uses a composition of smooth approximations towards the max- 
imisation of the data fitting, ensuring a good matching to the input data as well 
as optimal interpolation characteristics. The pipeline can choose between degree of 
data fitting and smoothness of the final meshes. The meshes used in this study were 
produced following the pipeline, presented in Figure A.1 and briefly described below, 
were obtained using a high degree of Laplacian smoothing in order to mimic heart 
muscle texture. 

The LV epicardium and endocardium, as well as the RV endocardium, were semi- 
automatically segmented using a digital drawing board and pen. The RV epicardium 
could not be resolved at this resolution and a uniform wall thickness of 3 mm of 
the RV wall was assigned, based on a large ex-vivo study of 30 ARVC hearts [167]. 
This is because the RV wall epicardial surface could not be distinguished from the 
pericardium in the MRI scans. As the scanner moves towards the apex of the ventricle, 
the uncertainty in the location of the RV epicardial surface grows as can be seen in 
Figure A.2. The mean RV wall thickness was reduced to 2 mm when investigating 
the effects of RV wall thinning. The short axis derived contours were aligned with 
the horizontal long axis and vertical long axis contours, by minimising the distance 


at expected intersection points in the plane of the short axis slices. 
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Figure A.1: Computational pipeline from MRI to personalised torso-biventricular 
meshes. A. Torso segmentation from MRI. The torso of an ACM patient was seg- 
mented semi-automatically (in red) by iso-intensity, selecting all points with the same 
intensity as the torso point selected, then removing any outlier points. B. Distribu- 
tion of torso contours for an ACM patient, colour-coded by acquisition view (coronal 
view in red, axial view in green, and short axis view in purple), and overlayed on two 
MRI images. C. Volumetric torso reconstruction. The original statistical torso model 
is shown (left), together with the maxima (+30’) and minima (—3c) of its first 3 
principal components, where o indicates standard deviation. The reconstructed torso 
is shown after the initial gradient descent fit, and after the subsequent thin plate 
splines deformation (middle and right panels, respectively). The distance between 
the mesh and contours was halved as a result. Contours are colour-coded by distance 
to the nearest mesh point (in millimetres). 


Three initial tubular meshes were then created for the LV endocardial, RV endo- 
cardial and whole epicardial surfaces. These initial meshes were fitted to the contours 
using a Thin Plate Splines algorithm [178] [179] which in brief applies small deforma- 


tions to the mesh and then measures the mean distance between each coarse mesh 
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Figure A.2: A-E. Increasing uncertainty in the location of the RV epicardial surface 
as the scanner moves from the base to the apex of the RV. The red arrow shows the 
location of the RV wall, visible in A-C in black and the RV wall and surrounding 
tissue merging in images shown in D-E. 


node and the closest contour point. The independent generation of the three surface 
meshes can lead to intersections of epicardial and endocardial surfaces in certain cases, 
especially in the thin RV. In that case, the intersecting surfaces are separated using 
an adaptation of the attractor-based mesh deformation algorithm. Attractor points 
were added near the intersection between the two surfaces. Epicardial surface points, 
which lie on the inside of the endocardial surface, become attractor points for the 
endocardial surface and viceversa. If this process, repeated twice over all intersecting 
points, was not sufficient to separate the two surfaces, then the endocardial surface 
still lying on the outer side of the epicardial wall became the epicardial surface and 
vice versa. Subsequent Laplacian smoothing, decimation and affine transformations 
ensured the new surface patches are well embedded in the surfaces. 

The space in between the surfaces is filled with tetrahedra once the three non- 
intersecting surfaces were generated. The goal is to achieve regular tetrahedra as this 
decreases numerical simulation errors [180]. An assessment of two of the most popular 
tetrahedra-making libraries, namely Gmsh [181]. The space in between the surfaces 
is filled with tetrahedra once the three non-intersecting surfaces were generated. The 
goal is to achieve regular tetrahedra as this decreases numerical simulation errors 


[180]. An assessment of two of the most popular tetrahedra-making libraries, namely 
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Gmsh [181] and Tetgen [182] was performed. The libraries performed similarly in 
most quality measures with the exception of the regularity of the tetrahedra created, 
measured by the ratio between the radii of the smallest possible sphere encompassing 
a tetrahedron and the largest sphere fitting entirely in it. Using this metric, Gmsh 
produced more regular tetrahedra and was used for volumetrisation. 

There are two major potential sources of uncertainty in the mesh creation pipeline. 
The first is the introduction of observer bias into the contours. The second source 
of uncertainty is the accuracy of fitting surfaces to a sparse set of contours. The 
second source of error was quantified using local uncertainty, defined as the maxi- 
mum distance between a contour and its nearest mesh surface area element. For our 
ACM cohort analysed, the mean contour point to mesh node discrepancy was 1.12 
mm across the 20 patients. This accuracy of the fitting is illustrated in Figure A.3 
for a representative ACM patient, evidencing minimal discrepancy between the con- 
tours of two intersecting acquisitions. This discrepancy was minimised by the surface 
generating algorithm, by choosing the middle ground between contours. 
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Figure A.3: Final volumetric mesh and initial contours for a representative ACM 
patient. Contours are colour-coded by distance to the nearest node in the mesh. All 
contours lie within 5 mm from the final mesh. The mesh is cut at the level of the 
most basal short axis slice contoured. 
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A.1.2 Construction of MRI-based torso meshes 


Images from the MRI cine stack, as well as the initial localiser slices, were selected 
and used in the torso contouring process. Slices were segmented in the areas of the 
shoulders, neck and the torso trunk beneath the arms. A subset of contoured images 
can be seen in Figure A.1A. Since MRI provides only a sparse set of images, volumetric 
reconstruction of the torso is challenging. Note the sparsity of the contours in Figure 
A.1B, which shows the complete contour set for a given patient. 

The proposed volumetric construction problem of the torso was solved using previ- 
ously developed methods [153] [183], by fitting a statistical male torso model (SMTM) 
generated from 4300 subjects [184] to the contours. After a rigid transformation, the 
first 40 principal component variations of the SMTM were used to fit the model to 
the contours. The initialisation of the torso position is in accordance with the posi- 
tion of the previously reconstructed heart from the patient. The optimisation used 
is a gradient descent-based algorithm, using the root mean square (RMS) distance 
between contour points and torso as a target function. Following this step, the result 
of the optimisation algorithm is further deformed by a thin plate spline algorithm, 
resulting in a halving of the RMS distance, as illustrated in Figure A.1C. 


A.1.3 Location of 12-lead ECG electrodes 


The electrode locations corresponding to the standard positions of the 12-lead ECG 
(used in the reconstruction of the QRS complex) were found by mapping the known 
electrode locations of a torso reconstruction onto the new reconstructed torso of each 
patient [161]. The mapping between the new patient torso and the torso with known 
electrode locations was found by means of a co-registration using thin plate splines 
deformation [185]. Using this mapping, the electrode positions were transferred onto 


the new torso of the patient. 


A.2 Simulation of the electrical activation sequence 


A.2.1 Ejikonal model of electrical propagation in tissue 


In the Eikonal equation, wavefront arrival times t, in the myocardium w are described 
based on a spatially heterogeneous orthotropic velocity function, V(x), and initial ac- 
tivations to at locations y, which are the connections of the His-Purkinje free running 
tree (not modelled explicitly here) with the endocardial ventricular surface in our 
model. The Eikonal equation is given by 
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Jt VVta = 1linw 
Where 
tg=t¢ in ¥ 


where w is the myocardial muscle and t, describes the wavefront arrival at location 
x. The symmetric positive definite velocity tensor V(x) holds the squared velocities 
in the transverse, longitudinal and normal fiber directions. These are the equivalent 
of the conductivities in the more complex bidomain or monodomain formulations of 
cardiac electrical propagation [160]. The points in y are the initial activation points 
which are activated at time to. Fibre orientations were assigned using a rule-based 
algorithm, varying from 60° to —60° in endocardial to epicardial layers with respect 
to the radial direction [24], based on canine left ventricle ex vivo measurements. In 
our control model of ventricular activation, anisotropic conduction was considered by 
assigning longitudinal, transverse and intersheet conduction velocities of 47.8, 31.5 
and 29.7 cm/s, respectively, in order to match the activation sequence dynamics of 
biophysically-detailed models of human ventricular conduction [153] [161]. 


A.2.2  Purkinje-endocardial root contact point localisation 


ECGI was used to inform the Purkinje-endocardial contact points in the RV of the 
patients, by matching their RV epicardial breakthroughs (Figure A.4). For ECGI, 
surface potentials were recorded in a previous study [139] at 256 electrode locations 
on the torso of the patients at resting heart rate, using a portable recording system 
(ActiveTwo, BioSemi, Netherlands), and used to reconstruct their epicardial break- 
throughs [139]. The electrodes closest to the standard 12-lead ECG positions, as 
localised on our torso reconstructions, were selected from the 256 electrode locations 
to reconstruct the patient clinical ECGs. A high/low pass Butterworth filtering, at 
0.5 and 35 Hz frequency cut-offs, were used to remove baseline wandering and high 


frequency noise, respectively. 


A.2.3. Pseudo-ECG calculation 


Pseudo-ECGs were derived from ventricular activations using the formula 


@.(x’) =D i (-VVn). Vand 
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Figure A.4: ECGlI-informed placement of Purkinje-endocardial contact points. 
Left: ECGI epicardial breakthroughs for the representative ACM patients of 
LV_LGE+/TAD-/R- and RV_LGE+/TAD+/R+ phenotypes, in both anterior and 
posterior views. Right: resulting epicardial breakthroughs in the control simulation 
for the same patients. 


LV 
RV 


where D is the diffusion coefficient of the torso, VV, is the spatial gradient of 
the transmembrane voltage, and r is the Euclidean distance between the electrode 
position x’ and a mesh node x. At each electrode position, the voltage generated by 


Te 


the ventricles was integrated at a given time and then repeated for all timesteps in 
the activation sequence, at intervals 1 millisecond apart. The diffusion coefficient D 
used was 0.1 m?/s for all the patients studied, as different torso conductivities on the 
ECG only results into a scaling of the amplitudes [161]. When a myocardial location 
was stimulated at time t,, its voltage evolution was approximated by a step function, 
with resting voltage of -90 mV and an activated voltage of 40 mV. 


A.3 LGE-MRI Informed Modelling of Myocardial Fi- 
brosis 


LGE-MRI scans were analysed by clinical radiologists as previously published [139], 
using the Otsu semi-automated technique [186]. The course ECGI mesh were seg- 
mented into 21 AHA-style regions by the radiologists [165] and a regional value of 
fibrosis was specified for each, depending on their fraction of enhanced area. In the LV, 
pixels of intensity above a given threshold were identified as areas of fibrosis, as pre- 
viously described [139]. In the RV, as the penetration of the contrast agent is poorer 
and hence the resolution is lower, the quantification of fibrosis was done visually by 
two independent experienced radiologists, and for each region either no/some/high 
contrast was reported. To transfer this information onto the finer simulation meshes 
(Figure A.5), the coarse and fine meshes were co-registered by the use of anatomical 
landmarks (centre of the LV and the RV insertion points), using the Thin Spline 
Deformation algorithm [178]. 
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Figure A.5: LGE-MRI data driven models of spatial fibrosis. A. Example of LGE- 
MRI scan for the representative LV_LGE+/TAD-/R- patient, highlighting a region 
of fibrosis in the RV free wall. B. Quantification of the AHA-segmented regional level 
of fibrosis. The lateral free RV wall and the anterior RV are heavily fibrotic. C. 
Resulting fibrosis region registered in our patient mesh. 
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